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Appendix A — LSF Cone Calorimeter Data 


A.1 Nomenclature 

A.2 Cone Calorimeter Data 
A.3 Ignition Data 

A.4 Heat of Combustion Data 
A.5 Heat of Gasification Data 
A.6 Total Energy per Unit Area 
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A. 1 — Nomenclature 


Due-to limitations in the presentation of the Cone Calorimeter data used, the 
following nomenclature shall apply to all tables and charts contained in this Appendix. 
All sample materials were tested 20 times with the cone calorimeter--exposure of 5 
samples to 25, 35, 40 and 50 kW/m’ external radiant heat fluxes. However not all 
samples ignited and only the values for materials that did ignite are presented in the 
following tables and graphs. The three average heat of combustion values represent a 
numerical average of the values obtained for specimens that ignited during testing. The 
following “symbols” are used to represent the actual physical values and have the units 
listed: 


Table A. 1: Nomenclature and Units Used for the Symbols Appearing in Appendices 
Graphs and Tables. 


Tet | a Wat td cae 
eG es 


The peak energy release rate of the material per 


Pan 4 Pan 4 
peak square meter. 


Opes 
Hc The heat of combustion associated with the peak | 
he ae ne energy release rate. 
i ee The specimen mass loss rate per square meter. 
M. "peak Lees g/s-m? an 
Calculated from O" 12.44 / AHC jong 


The average energy release rate around the peak 
Oar kW/m? 


release rate. The integrated average of the heat release 
rate (Q") values that are above 80% of O" 1.24 - 
The integrated average of the heat of combustion 
HCpeak avg. AH ~, peak avg. kJ/g 
calculated. 
3 The average specimen mass loss rate per square meter. 
M.” peak avg. MM sneak avg. g/s*m Calculated from Q" / AHc 
peak avg. peak avg. 
The negative value of the slope through the specimen 
-dm/dt —dm/dt g/s mass vs. time curve based on the region of actual 
burning. 
ey. or pe The overall energy release rate of the material 
- overall avg. age kW/m : on" 
overall avg per square meter: AH. Sane Spat pe 
The overall heat of combustion for the specimen. 
HCoverail avg. AH Coverall ave. kJ/g Based on the total energy evolved divided by the 
total mass loss. 
The overall specimen mass loss rate per square meter. 
=f 92 
Tt earch avg. g/ s°m 


values that are over the same range that O" is 
peak avg. 


Calculated from (- dm/ dt)/ A epic » Where Asampte = 
0.0088 m7. 


M." overall avg. 
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A.2 — Cone Calorimeter Data 
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Cone Calorimeter Data R 4.21 F.R. Expanded Polystyrene Board (80 mm) 
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Cone Calorimeter Data R 4.21 F.R. Expanded Polystyrene Board (80 mm) 
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Cone Calorimeter Data R 4.21 F.R. Expanded Polystyrene Board (80 mm) 
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Cone Calorimeter Data R 4.21 F.R. Expanded Polystyrene Board (80 mm) 
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Cone Calorimeter Data R 4.21 F.R. Expanded Polystyrene Board (80 mm) 
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Cone Calorimeter Data R 4.21 F.R. Expanded Polystyrene Board (80 mm) 
40 kW/m’, Test #1 
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Cone Calorimeter Data R 4.21 F.R. Expanded Polystyrene Board (80 mm) 
40 kW/m”, Test #3 


Heat Release Rate 


0 50 100 150 200 250 300 350 400 
Time (sec) 


Heat of Combustion 


Heat of Combustion (kJ/g) 


0 50 100 150 200 250 300 350 400 
Time (sec) 


Specimen Mass 


0 50 100 150 200 250 300 350 400 
Time (sec) 


Heat Release Rate (kVWV/ 


Heat of Combustion (kJ/g) 


Cone Calorimeter Data R 4.21 F.R. Expanded Polystyrene Board (80 mm) 
40 kW/m”, Test #4 


Heat Release Rate 


0 50 100 150 200 250 300 350 400 
Time (sec) ; 


Heat of Combustion 


0 50 100 150 200 250 300 350 400 
Time (sec) 


Specimen Mass 


00 150 200 250 300 350 400 
Time (sec) 


Cone Calorimeter Data R 4.21 F.R. Expanded Polystyrene Board (80 mm) 
40 kW/m’, Test #5 
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Cone Calorimeter Data R 4.21 F.R. Expanded Polystyrene Board (80 mm) 
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Cone Calorimeter Data R 4.21 F.R. Expanded Polystyrene Board (80 mm) 
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A.3 — Ignition Data 
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A.4 — Heat of Combustion Data 
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Downward and Lateral Flame Spread in 


Roland Apparatus Phase 5 


Dr. James G. Quintiere and Chen-Hsiang Su 


~ Abstract 


The test data developed from the Roland apparatus phase 5 are analyzed by using a 
theoretical flame spread model. The Roland apparatus considers a vertically oriented 
specimen with a 40 kW/m’ radiant heat source and a 0.42 kW pilot flame. The specimens 
were divided into 10cm by 10cm areas in order to observe the flame spread on the 
specimen. During the test the lateral and downward flame developments to these small 
areas are recorded. By plotting the lateral and downward flame spread velocities versus 
irradiance, an empirical correlation based on the theoretical model is found to correlate the 
test data. As a consequence, the material flame heating parameter, ® , is derived for further 


application. 


NOTATION 


parameter for thermal response function 
specific heat 

negative value of the slope fitting the test data 
thermal response function 

surface heat loss coefficient 

convective heat transfer coefficient 

thermal conductivity 


thermal inertia of material 


external radiant heat flux per unit time per unit area 
critical heat flux for flame spread 


critical heat flux for ignition 
time 

characteristic equilibrium time 
temperature 

ignition temperature 

surface temperature 

ambient temperature 

flame front velocity 


distance along the test specimen 


density 
flame heating parameter 


Stefan-Boltzmann constant 


1. INTRODUCTION 


In recent years, fire modeling technique has been recognized as a powerful and 
convenient tool for predicting the real fire phenomena. However, the accuracy and 
applicability of the fire model strongly depend on the sufficient understanding of related 
material properties. In this paper, a method of deriving material flame spread parameter ® 
from existing test data is proposed. The data used in this paper is offered by the LS Fire 
Laboratory (Italy) for typical construction materials which contains flame spread data from 


Roland Apparatus Phase 4 and 5. 


The Roland test [1] apparatus is designed to simulate a starting fire under condition 1.) 
Fire in a corner, bf condition 2.) Fire exposure on the wall/specimen of 40 kW/m’. The 
fire source is the ISO radiant panel which is positioned in an 35 degree angle to the 
specimen. The set up gives a similar feed back to the specimen as observed in a real corner 
fire. The time of ignition and the times which the lateral and downward flame fronts 


reached each 100mm increment were recorded and mapped. 


Under the thermally thick assumption, the theoretical flame spread velocity could be 
written as V = ® (kpc)” (Ti,-T,)” {2]. By applying the lateral and downward flame 
spread data from the Roland test, the parameters , kpc, and Ti, are deduced by 
analysis. Also a critical surface temperature is determined which continuous spread is not 
possible. It should be mentioned that the theoretical model does not consider melting, 
charring, regression, or the fluctuation of Ti,, therefor the deduced parameters should be 
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considered effective only when the realistic flame spread conditions are in well compliance 
with the assumptions. Inaccuracy of the analysis could be expected for materials with 
excessive melting and dripping characteristics. However, the analyzing process described in 


this study provides a consistent approach of relating the test data to theoretical application. 


2. THEORETICAL DEVELOPMENT 


The basis of the theoretical model is displayed in Fig-1. For the semi-infinite solid 


under external radiant heat flux, the surface temperature T, can be given [2] as 


Tee T= a. F(t) , (1) 
2h vt_ , t > small 
Vazkpc 
where F(t) = (2) 
] , t > large 


In the above equations, q , _ is the external radiant heat flux, T., is the ambient and 


initial temperature, h is the surface heat loss coefficient, and F(t) is a function of 


time and the thermal properties of the solid. 


1h Thick solid ae. 


Fig-1 Flame spread model 


For the steady state condition after long heating time, F(t)—>1, equation (1) can be 


rewritten as 


q eet Deciel sh G-Cl ates Ta) wee OTe) (3) 


where h, is the convective heat transfer coefficient. 


By substituting T,= Tj, into equation (3), the critical heat flux for ignition q , ig 


could be found as 


ap =he( Ty --T.)+0( To 1.4) = Ge ee (4) 


Based on the equation (1) to (4), an empirical result has been found to describe the 


ignition data. It can be written as 


ae -Fo=| popes ©) 
q. : 


V4 
+ 


> 


where b is aconstant fora given material, and t* isthe characteristic time indicating 
the thermal equilibrium for F(t*) = 1. By comparing equation (2) and (5), the following 
relation could be found. 


 Inkpc 


Equation (6) could be used to determine F(t) while h and kpc for a certain 


material are known. 


The flame spread velocity V_ is expressed [2] as 


® 
AF eee 7 
koe Gre ry)? i 


where @ is the flame heating parameter for a given material. Equation (7) can also. be 


written as 


Viati2 Cais é i“ (8) 


i 


By substituting equation (3) and (4) into equation (8), the equation becomes 


Vie ca ae e 4) [a 7 : Fo) 
k Pig an ER 
ales 3 sé cre) 
(DQ Bei) 


sila ie a: al (9) 


where @® = (4/n) (Cb) ' (10) 


and -C_ is the slope observed from the test data plotof V~? versus q , F(t). 


3. THE ROLAND TEST PHASE 5 APPARATUS 


3.1 DESCRIPTION 
The Roland test apparatus is designed to simulate a starting corner fire as shown in 


Fig-2. 


Real fire | 


Fig -2 Concept of Roland Apparatus 


The radiation panel gives a constant output of 40 kW/m which together with the 0.42 kW 
pilot flame create a stable environment for testing. The specimen to test is 1m by 1.5 m and 
is placed in a specimen holder which again is placed on a electrical driven trolley. The 
moving of the test specimen is controlled from outside the test room so as to ensure stable 


test conditions. 


The test specimen is placed on top of an insulated steel tray to capture the eventually 
dripping from the burning specimen. Capturing of dripping is important for the 


development of the fire and is therefor visually observed and recorded. 
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~WBWBWBALRARAABBRAABE 


Test storts when ddtaT > 10 C 


Fig-3 Test Apparatus Setup 


in front of the radiation panel to ignite the 


A movable pilot flame is mounted 


combustible gasses from the test products. The length of the flame is approximately 25 mm 


and they are not touching the surface of the product as shown in Fig - 4. If the product 
melts or expands, the pilot flame could move with the surface of the product so as to 


uJ 
: 
: 
1 
a 
g 
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Fig-4 The Pilot Flame 


simulate the real fire condition. 


3.2 TEST CONDITIONS 

The room for Roland apparatus is 19 m’ in area and has a volume of 74 m*. An 
extraction hood below the ceiling covering the specimen and the radiant panel is mounted 
in the room. The air change rate is 14.5 times/hr. 2 symmetrically placed openings at the 
end of the room serve the air and the air velocity was kept under 0.1 m/s so as to avoid any 


turbulence which might affect the burning behavior of the test specimen. 


3.3 INSTRUMENTATION 
The design of the extraction duct, the measuring device and the analyzers follow the 
specifications from the ISO 9705 Room Corner. The result was logged into a computer 


each 6 seconds. 


3.4 CALIBRATION OF THE RADIANT PANEL 
In order to secure the agreed level of the thermal attack of 40 kW/m’ in an area of 300 
cm on the specimen, the heat flux distribution was checked frequently and the result was 


shown in Fig - 5. 


3.5 TEST SPECIMEN 
The test specimen was prepared with a size of 1000 by 1500 mm, as in Fig - 6. The 
specimen was marked and divided into 100 by 100 mm areas. For a more detail dimension 


of the setup, see Appendix . 
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Fig-5 Roland Phase 5 Heat Flux Calibration Diagram 


Fig-6 Dimension of the test specimen 
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Fig-7 Illustration of the test apparatus 
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4. ANALYSIS PROCEDURE 


4.1 FLAME SPREAD TEST 


During the testing period, the measurements and observations include: 
Maximum temperature rise in the hood, 

Production of the heat release, 

Light intensity, 

Time to ignition, 

Time for flame front to reach 100 mm increment in upper, lateral, and downward 
direction, 

Duration of flaming, 

Critical flux, 


Dipping, burned area, damaged area, and weight loss. 


A sample flame spread record of the test result is shown in Fig - 8. 


PUR WITH PAPER GLUED ON CALCIUM SILICATE 5.04 


i 

| 
pepouras 

s @anZ 


Tome to igniten 
Grretean of Gaming 


. 
8 g = 


Ouration of teat 


> Average fase sprue epeexio quitin. 000 
© TOTALLY DESTROYED AREA Average Game spread Gowmwarts  awwimin. 3101.05 


‘S) OMMACLD ARLA (FLAMING ) Avengo temo sprusd teens euvemin, 974.16 


Fig-8 A-sample record of the test result 
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4.2 DATA ANALYSIS 

(1) Obtain kpc and Tj, from Roland apparatus phase 4 
From a separate study on the test result of Roland apparatus phase 4 test data [3], the 
kpc and Tj, for several materials were derived. Those derived values were used 


as input parameters for this study. 


(2) Compute flame front velocity V 
By applying a running three-point least square fit , as in equation (11), to the measured 
flame front position-time (x,t) data, the flame spread velocity could be found. 


se mt 2x 
oo (11) 
Sie (Xt) 
3 


(3) Calculate the surface heat loss coefficient h and critical heat flux q ,;, 
Since Tj, is known, equation (4) would than gives the value of q ,;, . The value 


of h, is chosen to be 15 W/m’, which is similar to the LIFT test [4]. After q ,;, 


was found, h_ is then derived by applying the equation (4). 


(4) Calculate the parameter b and F(t) 


Withthe kpc and Tj, previously derived from Roland apparatus phase 4, the 


parameter b_ could be found by applying equation (6). For any given t_ , F(t) 
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is then determined by using equation (5). 


(5) Find the slope -C 


As shown in Fig - 9, by plotting V'* versus q | F(t),theslope -C canbe 


estimated either mathematically or manually. The fitting line was forced to pass trough 


the point q ,;, sos to obtain a more consistent correlation. 


(6) Compute the flame spread parameter ® 
After the value C was determined in the previous step, the flame spread parameter 


@ is finally determined by using equation (10). 


yi2 


a | F(t) 


Fig-9 Estimation of the slope -C from derived data plot 
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5. RESULTS 


Based on the proposed analysis procedure, totally 13 materials were analyzed and a 
summary of the results is sown in Table - 1. Detail calculating for the materials that flame 
spread were observed are shown in Table - 2 to Table - 19. The tables show the results for 
varied external irradiances for no pre-heat conditions where the solid line indicate the slope. 
It should be noticed that the high end of the solid line is near extinction where chemical 
kinetic dominates and heat loss phenomena is suspected to influence the flame spread speed. 
The low end, on the other hand, is near ignition where spread velocities are high and 
ignition phenomena can effect flame spread. Therefor, scattered data points at the both 


ends of the solid line could be expected. 


For the wooden material group, 5-04, -09, -10, -11, the test data correlate well with 
the fitting lines, and the magnitude of derived flame heating parameters are within 
reasonable range. For the thermal plastic group, 5-05, -07, -08, -12,-13, test data were 
dispersed and did not fit well with a straight line. There were no flame spread at all for 
three materials, 5-01, -02,-03. There is also no test data for material 5-06 due’ to melting 


and unusual flame spread. 


While the properties of thermal plastic materials do not tend to follow the theoretical 
assumptions as described in the introduction, it is not surprised that the data were dispersed. 
There are two possible major effects that caused the difficulties of correlating the data, one 


is the melting effect, another is the tray effect. When heating up the products in this group, 
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the flame spread velocities are depending on the melting-away patterns, dripping levels, 
and the continuation of the burning of the droplets after they have been collected in the tray. 
For cases that the data were obviously influenced by both of the two effects, two ®s 


were derived for reference. 


It is unusual that the flame spread downwards but not laterally for material 5-08. This 
could be caused by the distortion of the specimen that blocked the flame propagation 
during heating up. In the general, since there are more measuring points at the lateral 
direction, the resolutions are better than the downward, and thus the estimation of @ is 
more accurate. For materials that only few test data are available, it is difficult to achieve 
accurate analysis results. Further investigation and experiments are required to reveal the 


related properties of those products. 
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6. CONCLUSIONS 


The analysis and results of this study illustrate a well defined procedure of deriving 
material flame spread properties from Roland test data. The flame heatmg parameters 
found in this study can be used either in advanced fire models or as the performance 


indexes of material fire resistance. 


A spread sheet program for analyzing the Roland test data is successfully built. With 
proper input data, the spread sheet will execute all require calculations and generate related 
data plots. Since the diagrams are directly linked to the data cells in the spread sheet, the 


estimation of the slopes could be done in few trial-error steps. 


Although only the test data of Roland apparatus phase 5 are analyzed in this study, the 
concept of proposed method is not limited to the specific apparatus. With a little 
modification of the analysis procedure, the method or the spread sheet program could be 


applied to similar test apparatus and give satisfied results. 
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Table - 1 Summary of Analysis Results 


Tig kpc  |@Lateral |{® Downward| Ts, min 
e (°C ) Ickw2m*eK?S)} (kW?2/m?) (kW ?/m?) tia ee) 


5-05 XPS 40mm l. 2 
oP a a PD 


S05 Lend ee 205 oo 


pis aes 
J 200.5 ** 


Notes: 
"+": No flame spread was observed. 
++": Not able to test due to melting and unusual flame spread. 
"*" + Tray effect suspected. 
"*#" = Not accurate due to scattered test data. 
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PUR WITH PAPER GLUED ON CALCIUM SILICATE 5.04 : 


2™°SPECIMEN 


100 80 60 70 00 6 40 W 20 10 0 


(48) TOTALLY DESTROYED AREA 


©) DAMAGED AREA (FLAMING ) 


SUMMARIZING LIST 


Orpping ( lev. ) 1 
Conditioned weight (or.)  12235.5 


Weight after test (gr. ) 11624.1 
(gr. ) 611.4 
(%) 5 
25.5 


Maximum flame spread upwards 
Maximum flame spread downwards 
Maximum flame spread left 


Maximum flame spread right 


Critical flux upwards 
Critical flux downwards 


Critical flux left 


Time to Ignition 
Duratlon of flaming 


Duration of test 


Average flame spread upwards 


Average flame spread downwards 


Average flame spread latoral 


mm. 
mm. 
mm. 


mm, 
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mmmin. 


mmmin. 
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4™SPECIMEN 


70 60 80 40 30 20 10 0 
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@) TOTALLY DESTROYED AREA 


(\) DAMAGED AREA ( FLAMING ) 


SUMMARIZING LIST 


Dripping (lev. ) 1 
Conditioned weight (gr. ) 12781.2 


Weight after test (or. ) 12234.3 
(gr. ) 546.9 
(%) 43 
23 
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Maximum flame spread upwards 
Maximum flame spread downwards 
Maximum flame spread left 


Maximum flame spread right 


Critical flux upwards 
Critical flux downwards 


Critical flux left 


Time to ignition 
Duration of flaming 
Duration of test 


Average flame spread upwards 


Average flame spread downwards 


Average flame spread lateral 


mm/min. 


mmmin. 


mnavmin. 
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PUR WITH PAPER GLUED 


6™SPECIMEN 


CB) TOTALLY DESTROYED AREA 
()) DAMAGED AREA ( FLAMING ) 
SUMMARIZING LIST 


,; i ( lev.) 
(gr. ) 


(gr. ) 
(gr. ) 


ON CALCIUM SIL 


nee lies! eh sn  E.., 


ICATE 5.04 


Maximum flame spread upwards 


Maximum flame spread downwards 


Maximum flame spread left 


Maximum flame spread right 


Critical flux upwards 
Critical flux downwards 


Critical flux left 


Time to ignition 
Duration of flaming 


Duration of test 


Average flame spread upwards 
Average flame spread downwards 


Average flame spread lateral 


_——— iL, 


mm. 1000 
mm. 500 
mm 700 
mm 100 
kw/m2 <2.60 
kw/m2 <9.70 
kw/m2 <2.10 
s 2 
$s 723 
8 1321 
mmmin. §550 


mimin. 1781.33 


mavmin. 328.23 
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Appendix C - Thureson, Per, “Fire Tests of Linings According to 
Room/Corner Test, ISO 9705”, Swedish National 
Testing and Research Institute, Fire Technology, Report 


95R22049, January, 1996. 
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Issued by an Accredited Laboratory 


L.S.F Fire LABORATORIES 
Via Garibaldi, 28/A 


22070 Montano Lucino (CO) 

ITALY 
Handlaggare, enhet/Handled by, department Darum/Date Beteckning/Reference Sida/Page 
Per Thureson, Fire Technology 1996-01-26 95R22049 1 (5) 


Tel +46 (0)33 16 50 83 


Heat release and smoke production of building products 
when tested according to the Room/Corner Test, 
ISO 9705 


(14 enclosures) 


SP - Fire Technology has by the order of LSF performed a series of tests of 
linings according to the Room/Corner Test, ISO 9705. The purpose of the tests 
was to evaluate the burning behaviour of the products. 


Products 
Twelve building products were tested. They are described in enclosure 1. 


The products were delivered to SP during November-1995. Products no 1-8 
were delivered from DBi, Denmark. Products no 9-12 were delivered by LSF. It 
is not known to SP if the products received are representative of the mean 
production characteristics. 


Test procedure 


The Room/Cormer Test - ISO 9705: 


The test room is constructed from aerated concrete and has the following 
nominal inner dimensions; 3.6 m x 2.4 m x 2.4 m (length x width x height). 
There is a doorway in the centre of one of the 2.4 m x 2.4 m walls. The material 
is mounted on the walls and in the ceiling of the test room. A gas burner having 
the dimensions 170 mm x 170 mm x 145 mm is placed on the floor in the right 
comer opposite the doorway, see figure 1. 


The smoke produced by the burning specimen is collected by the exhaust hood 
and evacuated through the exhaust duct. By continuously analysing the oxygen 
(O2) and carbon dioxide (CO?) content of the smoke gases, the developed heat 
release rate can be calculated. The optical density of the smoke is measured 
continuously throughout the test with a white light system. 


The test duration is 20 minutes. During the first 10 minutes the heat output from 
the burner is 100 kW. Thereafter the output level is increased to 300 kW for 
another 10 minutes. If flashover occurs the fire is extinguished. 


SP, Sveriges Provnings- och Forskningsinstitut, Box 857, 501 15 BORAS, Tel 033-16 50 00, Telefax 033-13 55 02, E-mail info@sp.se, Org. nr 556464-6874 
SP, Swedish National Testing and Research Insttute, Box 857, S-501 15 BORAS, SWEDEN, Telephone + 46 33 16 50 00, Telefax + 46 33 13 55 02, E-mail info@sp.se, Reg.No 556464-6874 


Laboratorium ackrediteras av Styrelsen for ackreditering och teknisk kontroll (SWEDAC) enligt svensk lag. Verksamheten vid de svenska ackrediterade laboratonema uppfylier kraven i 
SS-EN 45001 (1989), SS-EN 45002 (1989) och ISOMEC Guide 25 (1990:E). Denna rapport far endast aterges i sin helhet, om inte SWEDAC och SP i forvag skrittligen godkant annat. 


Laboratories are accredited by the Swedish Board for Accreditation and Contormity Assessment (SWEDAC) under the terms of Swedish legislation. The Swedish accredited laboratories meet the 
requirements in SS-EN 45001 (1989), SS-EN 45002 (1989) and ISONEC Guide 25 (1990:E). This repon may not be reproduced other than in full, except with the prior written approval of 
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Optical density ——Gas analysis (O02, CO, CO) 
(lamp/photocell) Ry stame flow 
| Temperature and differential pressure 
Exhaust hood 


3.0m x 3.0m x 1.0m 


2.40m 


Doorway 0.8m x 2.0m 


Figure 1 Schematic drawing of the ISO 9705 room test. 


Surface temperatures: 


The surface temperatures of the test samples in the ceiling and on the right wall 
(viewed from the door opening) were measured with thermocouples, 0.25 mm 
thick lead wires of type K (NiCr-NiAl). Five thermocouples were mounted in 
the ceiling and another five were mounted on the upper part of the wall. The 
positioning of the thermocouples was according to figure 2. The hot junction of 
the thermocouple was fixed in contact with the surface of the wall/ceiling by 
means of non combustible tape. The lead wires were led through the wall and 
ceiling of the room, close to the hot junction. 


Figure 2 Positioning of thermocouples in the ceiling (left) and on the wall 
(right) (measures in mm). 
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Heat flux measurements: 


Heat flux measurements were made using total heat flux gauges that were facing 
the burner corner. Two heat flux gauges were mounted along the diagonal of the 
room at a distance of 1m and 2 m measured from the burner corner. The gauge 
at 2 m distance was mounted at 1.2 m above floor level. The other gauge, closer 
to the burner corner, was mounted at 0.8 m above floor level. The sensing 
surfaces were vertically oriented and facing directly the burner corner. 


Test results 


Tabulated test results are given in enclosure 1. Detailed information is given for 
each test in enclosures 2-13, including 


-logs of test observations 

-graphs of heat release rate (HRR) 

-graphs of smoke production rate (SPR) 
-graphs of carbon monoxide production rate 
-graphs of surface temperatures 

-graphs of heat flux 

-photographs 


Test results from a calibration test are given in enclosure 14. 


Summary and comments 


The performed test series showed that ISO 9705 is suitable for ranking products 
in a scientific sound way. Measured heat release rates from all the products are 
shown in figures 3-4. It is seen that the method accurately can measure the 
burning behaviour of very different products. The hood system used was 
capable of very accurate measurements ranging from 5 kW to 2000 kW (see also 
enclosure 14). Thus even very limited combustion can be detected. 


: 


Heat release rate [kW] 


Time [min] 


Figure 3 Products that reached HRR peaks exceeding 1000 kW . 
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Figure 4 Products that showed limited burning i.e. had HRR peaks lower than 
1000 kW. 


The ignition regime of the ISO 9705: 


The need of the second, more severe, heat output level (300 kW) of the burner to 
achieve the required discrimination is clearly demonstrated by comparing the 
heat release rate history of for example FR plywood and Gypsum board, see 
figure 5. Until the increase in burner heat output the behaviour of the products 
are quite similar, while the FR plywood shows a rapid fire growth rate after the 
burner heat output is increased. 


ses 245 


Heat release rate (kW] 


3 


Time [min] 


Figure 5 A comparison of the heat release rate history of FR plywood and 
Gypsum board. The 100 kW burner level is hardly capable of 
separating the burning behaviour of those products. However after 
the increase in burner heat output to 300 kW the difference in 
combustion properties is clearly seen. 
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The same conclusion can be drawn from the three polystyrenes test results. 
During the first ten minutes of the tests the polystyrenes behaved quite similar 
with HRR peaks at about 1000 kW, see figure 6. 


Heat release rate [kW] 


Orr 5 10 15 20 fh) 
Time [min} 


Figure 6 Heat release rates of the three tested polystyrenes. The burning rates 
are quite similar during the first period of the tests. However after 
increase of the burner heat output the differences in burning 
behaviour are clearly seen. 


However after the burner heat output was increased at 10 minutes two of the 
polystyrenes showed more than 1000 kW HRR but clearly separated in time. 
The third product showed only limited burning. 


Sveriges Provnings- och Forskningsinstitut 


Fire Technology, Materials Reaction-to-Fire 
Seti 


Bj¢@m Sundstr6m Per Thureson 
Head of Section Technical officer 


Enclosures 


1 Tabulated test results and list of tested products 

2 Paper-faced gypsum board, test results 

3 FR PVC, test results 

4 Acrylic glazing, test results 

5 FR extruded polystyrene board, 40 mm, test results 
6 PUR foam panel with aluminium faced paper finish, test results 
7 Mass timber, varnished, test results 

8 FR chipboard, test results 

9 3-layered FR polycarbonate panel, test results 

10 FR expanded polystyrene board, 40 mm, test results 
11 FR expanded polystyrene board, 80 mm, test results 
12 Plywood, test results 

13 FR plywood, test results 

14 Calibration test results 
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Enclosure 1 


_ Products and test results, tabulated data 


Table 1 Tested products 


p[Prsteeee [Res tm | eas) 
thickness [mm] | [kg/m3] 
[6 [Mass timber (ine), varnished [as 
[a | Sayered Fk polyeatbonat panel [16 [1200 
[9 FR expanded polysyrene board [030 
[10 [FR expanded polysrrene board [80 [15 
[Tiny [Rios ag ais eae aia ella ly 
[1a [ FR plywood kr aman an [ara ssp mama oe 


*No data given by the client. Measured data are given in enclosures. 


Table 2 Heat, smoke and carbon monoxide production. 


Time to 1000} Peak HRR Peak SPR 
Type of material kW or time to] excluding {m2/s] 
peak HRR if | the burner 
the fire was output 
smaller than [kW] 
1000 kW 
{min:s] 


Peak CO 
[g/s] 
Gypsum board 
FR PVG SSL Ts a ce 
Kenic glaze |e [tO 8 


FR ox, PS 40 mim 


PUR foam withalpap |__| >100_[ Sa 


Nias dimer a 
20:00 


11:36 
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Enclosure 1 


Table 3 Heat, smoke and carbon monoxide production, accumulated 
data (all data are calculated from start of test until the HRR reached 
1000 kW or to the end of the test duration) 


TSP Tot CO 
Type of material [m2] [g] 


Gypsum board si board 


RRC eae 
[Ses a oA FS A 
FRew Psamm | 7 | ™ | ™s [ 2 | 
PURfoamwihalrep | 2 | 5 | © | "| 
FR ehipbowrd |__| eR 
Bayer PR povetbon [WO [nos | zee 


TSP/THR 
[m2/MJ] 


The polystyrenes behaved similar during the first part of the test. A peak HRR just 
exceeding 1000 kW was seen quite early, see figure 6 (main report). However 
after the first peak HRR the fire decreased for all polystyrenes and a second peak 
occurred later on after the increase of the burner heat output. Therefore data are 
also given for the polystyrenes for the entire test duration, see table 4. 


Table 4 Polystyrenes, heat, smoke and carbon monoxide 
production. The values are given until the second time the HRR 
reached 1000 kW or to the end of the test duration (EPS 40 mm) 


the test or 
to end of test 


(Mj) 


ESSE) CH 
5 ER 
Epsom me eos [ese 
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Peak SPR: 


TSP: 


Peak CO: 


Tot CO: 


TSP/THR: 
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Enclosure 1 


Time at which the total heat release rate reached 1000 kW (including 
the burner heat output) corresponding with flashover in the room. 


Maximum heat release rate for the entire test duration if no flashover 
occurred (excluding the HRR from the ignition source). For tests 
where flashover occurred this value cannot be given since the fire 
had to be extinguished. 


Total heat released during the entire test for tests where no flashover 
occurred (excluding the HRR from the ignition source). For tests 
where flashover occurred the heat release rate is integrated from the 
beginning of the test till flashover. 


Maximum smoke production rate for the entire test duration if no 
flashover occurred. For tests where flashover occurred this value is 
given at flashover. The smoke production rate is calculated as 
follows: 


SPR = me Lee Vo 
L I 
Where 


V is the volumetric flow rate in the duct at actual temperature (m3/s) 
L is the optical path length in the duct (m) 

To is the initial intensity of a light beam 

I is the intensity of the light beam after traversing a smoky envi- 
ronment 


Total smoke production during the entire test for tests where no 
flashover occurred. The total smoke production is integrated from 
the smoke production rate. For tests where flashover occurred the 
smoke production is integrated from the beginning of the test till 
flashover. 


Maximum CO production rate for the entire test duration if no 
flashover occurred. For tests where flashover occurred this value is 
given till flashover (including any contribution from the ignition 
source). 


Total CO production during the entire test for tests where no flash- 
over occurred. The total CO production is integrated from the CO 
production rate. For tests where flashover occurred the CO prod- 
uction is integrated from the beginning of the test till flashover. 


Total smoke production divided by the total heat released. 
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Table 5 Heat flux measurements. The heat flux gauges were 
mounted along the diagonal of the room, facing the burner corner, 
at a distance of 1 m and 2 m measured from the corner. 


Material 


Gypsum board i; 

FRPVC PE PE EEC 
Aewlic gaamg | FO] FO | 329 7s 
FRenmded PS@0mm | FO | FO | 3] 172 _| 
a ED EE 
FR chipboard «247 234_[ _NA_| NA 


Bayer Fe poijenonse | 186 [vas | wA [WA 
Remade mn [FOF [ss se 
og Sea (areae| Parma | 
paar ac er RRM he 


Legend 

FO: Flashover 

NA: Not appropriate. 

Ea No heat flux meters were installed. 

Maximum Maximum heat flux during the entire test recorded by a total heat 

Heat flux flux meter positioned 1 m from the burner (including the 

1m: contribution from the ignition source). Tests where flashover 
occurred are not reported. 

Maximum Maximum heat flux during the entire test recorded by a total heat 

Heat flux flux meter positioned 2 m from the burner (including the 

2m: contribution from the ignition source). Tests where flashover 
occurred are not reported. 

Maximum Maximum heat flux till flashover recorded by a total heat flux 

Heat flux meter positioned 1 m from the burner (including the contribution 

Imtill FO: from the ignition source). In most cases this will be the heat flux at 
flashover. 

Maximum Maximum heat flux till flashover recorded by a total heat flux 

Heat flux meter positioned 2 m from the burner (including the contribution 

2m till FO: from the ignition source). In most cases this will be the heat flux at 


flashover. 


- Test results, ISO 9705 (NT FIRE 025) 


Product 
Paper-faced gypsum board. 
Mounting 


DanunvDate Beteckning/Reference 
1996-01-26 95R22049 


Sida/Page 
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Enclosure 2 


The gypsum boards were nailed to the lightweight concrete walls and ceiling. 


Observations during test. 

Time, [min:s] Observations 

0:00 Ignition of the burner, 100 kW. 

0:00-10:00 Very limited HRR and SPR. The paper surface was charred 
to an extent of approximately 1.5 m2. 

10:00 The burner output was increased to 300 kW. Some flame 


spread was seen on the ceiling after increase of the burner 


SP Sit 


heat output, see photo no 3. 


20:00 The test was terminated. The paper surface was charred 
to an extent of approximately 5 m2. 

After test: Charred areas, see photo no 5 and 6. 

Measured data 

Thickness, mm 12.6 

Density, kg/m? 720 

Conditioning 


Temperature 20+5 °C 
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Test results, graphs 


350 


Heat release rate [kW] 


0 5 10 15 20 Pie 30 
Time [min] 


Smoke production rate [m?/s] 


0 5 10 1S 20 25 30 
Time [min] 


CO production rate [g/s] 


0 5 10 Lo 20 23 30 
Time [min] 


Figure 1 Paper-faced gypsum board, heat release rate (including the HRR from the 
ignition source), smoke production rate and carbon monoxide production rate 
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Ceiling temp. [°C] 


0 5 10 15 20 25 30 
Time [min] 


Wall temp. [°C] 


es aa ER eles, Lobe Neal Mii Sass WA SUN Et ER 


= 
-— — 
. 


0 3 10 aS 20 Zs 30 
Time (min] 


Heat flux [kW/m?} 


0 5 10 15 20 25 30 
Time [min] 


Figure 2 Paper-faced gypsum board, ceiling temperatures, wall temperatures and 
heat flux, including the contribution from the ignition source. 
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Photo no 1 Prior to test "Paper faced gypsum board" ; 


The paper faced gypsum board were nailed to the lightweight concrete walls and 
ceiling. 


Photo no 2 Time 2:14 "Paper faced gypsum board" 


Very limited HRR and SPR (by request from the client standard polyether foam 
blocks were placed at three positions on the floor prior to start of test). 
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Photo no 3 Time 10:07 "Paper faced gypsum board" 


The burner output had been increased to 300 kW. Some flame spread was seen in the 
ceiling after increase of the burner heat output. 


"Paper faced gypsum board" 


Photo no 4 Lime 15732 
Stull very limited HRR and SPR. 
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Photo no 5 After test "Paper faced gypsum board" 


Most parts of the ceiling and the walls were undamaged. The paper surface was burnt 
or charred to an extent of approximately 5 m?. 


Photo no 6 After test “Paper faced gypsum board" 


Damaged areas in the ceiling. 
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Enclosure 3 


Test results, ISO 9705 (NT FIRE 025) 
Product 

FR PVC 

Mounting 


The PVC panels were screwed to a frame of light steel profiles which gave a 
spacing of 40 mm to the light weight concrete walls and ceiling. 


Observations during test. 


Time, [min:s] Observations 

0:00 Ignition of the burner, 100 kW. 

0:30 The ceiling panel above the burner began to deform, see 
photo no. 2. 

1:10 A smoke gas layer was formed. The ceiling was no longer 
visible, see photo no. 3 

J EF 2s) The material melted at the burner corner. 

9:00 Most of the ceiling material was melted (softened) and had 
fallen down onto the floor. HRR and SPR were limited 

10:00 The burner output was increased to 300 kW. 

11:00 All ceiling panels were softened and had fallen down. 

20:00 Gas shut off. 

After test: Not much material had been consumed. Most wall and 


ceiling panels had melted or just softened and were laying 
on the floor, see photo no 6. 


Comments to the graphs 


The accuracy of the HRR measurement was +10% due to noise pick-up in the 
analysers. 


Measured data 

Thickness, mm 3.0 
Density, kg/m 1505 
Conditioning 


Temperature 20+5 °C 


SP 511 
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Test results, graphs 


Heat release rate [kW] 


eocceccovecsersen}overeccseraccessoesces: 


Smoke production rate [m?/s] 
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Time [min] 


CO production rate [g/s] 


Figure 1 FR PVC, heat release rate (including the HRR from the ignition source), 
smoke production rate and carbon monoxide production rate 
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Figure 2 FR PVC, ceiling temperatures, wall temperatures and heat flux, 
including the contribution from the ignition source. 


,e2¢ ~e, ; ; 
‘ hen $> git tT \ rime ii ‘eel eit + ~ 


o-* 
2 * " } v3 wee , 
A Ss « 5 ig 
-— | .~ i ‘dei 
oe - * mg e Tiles pee ot mh ge a att eee sea De aed Wr» 
: be i* ie | q* r" a 
F FI pte os oe es 


: pueciona 


a y ‘ Py c-. a 
i a " - 
z ‘ 3! ' “a : 
: | 
® 1 Vy ' =" 
on 
. : } ' 
, @st results, gray hs phe ee eg ; ; 

op td wth ras Tes Tere Weer Pee Wes Oe Na 

\ 

f ; ~* 
— o=-otihepaniel ae ie “e er at veggie ‘ 

4 eo i be us ay 

e 24 e, ~. i | 
~—~—- — oe > “ aces eed 

; alae res 4 fend 
pH eS - m die «9 we aro ae 

‘. 7 

: (lice, ee ae ica, 

ae Say Ak idieniinete ccna uae <4 do laortatd aT 
fa 2 


a » of 
©. «Tee — oo + eee | gre epeseme aes i we — 


antes Sia oar ee rn = aye ee 


' ' 


— ™ a aw oe a a, a hem bet mp wh eet tenes + onde eee 


AfAm sp lain : ee ed 


— — - hon nail Tome " Sb ay : 
- ; Ae Seer tt an a | we a a ere ae 


~— —' _ 5 : spe ys es it ° ‘ . e how 


te ~~ ap ter ope epee ray Hanh Ove 
~ ee 2 . a 
& ’ a em ‘ ' i 
he 5 ae - ’ es . ’ | ein, Be ole Napmmee 9), nh f — bow seve 
j J, : nats ren ee 
/ +a ‘ ‘ ; 
} , ‘ ' : Pomel * , i 
-“—— a re ll eye ea ra J ct 4 Se eee ae oan 
4 ? ° ' . 


ai 


ik} : We 
1h eit oN zai _ ) Eh 
. eer 


* 
y 


i ae 
~ ¥ 
; va 
oh 
; * 
7 


2 owes ke Py. rer 

' ; fy oe 7 a 

r q* 5 : 5 

. a) \ 1 ; i ® " is - 
~- te ee some ate a — ~ 4 ~- i 

. : es he sehen par * 


ee ee ; 
a 
» ho yt tame ond ay tie ee rer 
‘© 2 


IP tae en 2 


zt . = ae a afl Bc 9 . ma) 
{aul sii IP Bia Ae ea 5: i eS “iat AS) 


. es Pare 

peennite : stow anna age 1 tN, 
) 4 eye a . ee ie eee 
Figure LER PVG Neat enienee re 


| amped Mats BEBADT 
,e¢ Spicnsw oe eas DatunvDate Beteckning/Reference Sida/Page 
se 1996-01-26 95 R22049 4 (6) 
Enclosure 3 


Photo no 1 Prior to test "FR PVC" 


The PVC panels were screwed to a frame of light steel profiles which provided for a 
spacing of 40 mm to the light weight concrete walls and ceiling. 


Photo no 2 Time 0:32 "FR PVC" 


The ceiling panel above the burner started to get deformed. 
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Photo no 3 Time 1:06 “EReEV i 


A smoke gas layer was formed. The ceiling was no longer visible. 


Photo no 4 Time 10:08 "FRPVC" 


The burner heat output had been increased to 300 kW. Large portions softened and 


fell down on the floor. Flaming in the material were only seen in the vicinity of the 
burner. 
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Photo no § Time 19:58 PERE Vic 


Close to gas shut off. Limited HRR and SPR, burning only in the burner corner. The 
material in the ceiling and most material on the walls had softened and fallen down 
on the floor. 


Photo no 6 After tes cPIcry GC 


Not much material had been consumed. The ceiling panels and most wall panels had 
melted or just softened and were laying on the floor. 
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Test results, ISO 9705 (NT FIRE 025) 
Product 

Acrylic glazing 

Mounting 


The acrylic glazing panels were screwed to a frame of light steel profiles which 
gave a spacing of 40 mm to the light weight concrete walls and ceiling. 


Observations during test. 

Time, [muin:s] Observations 

0:00 Ignition of the burner, 100 kW. 

1:00 ee in the burner comer were ignited, see photo 
no. 2. 

1:20 Melted material started to form a pool fire near the burner. 

1:30 The ceiling above the burner had ignited, see photo no. 3. 

1:45 Several burning droplets were formed. 

240 Approximately 2 m2 of melted acrylic was burning on the 
floor. Limited smoke production. 

2:39 Flames out the doorway. 

3:00 Extinguishment. 

After test: Most material was burnt or melted. The melted material 


covered approximately 50% of the floor. 
Comments to the graphs 
After two minutes the ceiling thermocuple "c1" (closest to the burner) was not 


anymore in contact with the ceiling material due to the melting of the material 
and hence "cl" measured only gas temperature from that point in time. 


Measured data 

Thickness, mm 3 
Density, kg/m3 1150 
Conditioning 


Temperature 20+5 °C 
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Heat release rate [kW] 
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Figure 1 Acrylic glazing, heat release rate (including the HRR from the ignition source), 
smoke production rate and carbon monoxide production rate 
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"Acrylic glazing" 


Photo no 1 Prior to test 


The acrylic glazing panels were screwed to a frame of light steel profiles which 
provided for a spacing of 40 mm to the light weight concrete walls and ceiling. 


"Acrylic glazing" 


Photo no 2 Time 0:57 


The panels in the burner corner had ignited. 
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Photo no 3 Time 1:38 “Acrylic giazing” 


The ceiling above the burner had ignited. A pool fire was formed. 


Photo no 4 Time 2:03 "Acrylic glazing" 


Several burning droplets were formed. Approximately 2 m? of melted acrylic was 
burning on the floor. Limited smoke production. Large portions of the walls and the 
ceiling had already been deformed. 
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Photo no 5 Time 2:27 "Acrylic glazing" 


Flash over was reached. HRR increased. SPR was limited. Burning ceiling material 
was falling down on the floor. 


Photo no 6 After test "Acrylic glazing" 


Large portions of melted acrylic had melted and run down on the floor. 
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Test results, ISO 9705 (NT FIRE 025) 

Product 

FR extruded polystyrene board, 40 mm 

Mounting 

The polystyrene boards were glued to a non combustible board called “Promatek 
H", density 870 kg/m3, with a water based contact adhesive called “Casco 


3880". The non combustible boards were nailed to the light weight concrete 
walls and ceiling before the polystyrene boards were glued. 


Observations during test. 

Time, [min:s] Observations 

0:00 Ignition of the burner, 100 kW. 

0:20 The material was ignited, see photo no. 2. 

0:50 SEP increase, the ceiling was no longer visible, see photo 
no. 3. 

1:25 HRR and SPR started to increase rapidly, see photo no. 4. 
The ceiling material melted. Burning droplets were formed. 

1:40 HRR peaked above 1000 kW, high SPR, see photo no. 5. 

2:00 HRR and SPR decreased. 

3:30 Only flames in the bummer corer. 

8:00 ' Melted polystyrene formed a pool fire near the burner. 

10:00 The burner output was increased to 300 kW. 

10:10 HRR and SPR increased rapidly. 

10:20 Flame spread down the walls, flames out the doorway. 
Burning droplets were formed, see photo no. 6. 

10:40 Gas shut off, extinguishment. 

After test: Most material was burned or melted, see photo no 8. 

Measured data 

Thickness, mm 40 

Density, kg/m? 33 

Conditioning 


Temperature 20+5 °C 
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Test results, graphs 
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Heat release rate [kW] 
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Smoke production rate [m/s] 
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Figure 1 FR extruded polystyrene board, heat release rate (including the HRR from the 
ignition source), smoke production rate and carbon monoxide production rate 
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Photo no 1 Prior to test "FR extruded polystyrene board, 40 mm" 


The polystyrene boards were glued to a non combustible board. The non combustible 
boards were nailed to the light weight concrete walls and ceiling. 


_ EXTRUDED 
POLYSTYRENE 


Photo no 2 Time 0:21 "FR extruded polystyrene board, 40 mm" 


The material had ignited. 
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Photo no 3 Time 0:57 "FR extruded polystyrene board, 40 mm" 


SPR increased, the ceiling was no longer visible. 


EXTRUDED 
POLYSTYRENE 


Photo no 4 Time 1:27 "FR extruded polystyrene board, 40 mm" 


HRR and SPR started to increase rapidly. The ceiling material melted. Burning 
droplets were formed. 
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_ EXTRUDED 
POLYSTYRENE 


Photo no 5 Time 1:42 "FR extruded polystyrene board, 40 mm" 
HRR peaked above 1000 kW, high SPR. 


Photo no 6 Time 10:23 "FR extruded polystyrene board, 40 mm" 


Flash over was reached. HRR and SPR increased rapidly. Flame spread downward 
the walls, flames out the doorway. Burning droplets were formed. 
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Photo no 7 Time: 10332 "FR extruded polystyrene board, 40 mm" 


HRR peaked above 2 MW. The entire room was engulfed in flames. Excessive SPR. 
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Photo no 8 After test "FR extruded polystyrene board, 40 mm" 


Most material was burnt or melted. The entire floor was covered with melted 
polystyrene. 


aah 
> Auacion 


Fy On 4 \ Dy OR 
u Be ph ay 


‘iu oa 
oF - 


pi Erk: ae . 
ue ae a 


sj heer At" 


79 big : ae , 


“mn GW brbed hieetniriok ake a" ) 


hot! oeh Aitw b boevey ah pi niles 


Ae 


Hi i . ii Cette 
aid : 


ran \ 
t 


. 


PR 
; ‘ 


SP 611 


=z Datum/Date Beteckning/Reference Sida/Page 
1996-01-26 95R22049 1 (6) 
Enclosure 6 


Test results, ISO 9705 (NT FIRE 025) 

Product 

PUR foam panel with al faced paper finish. 

Mounting 

The PUR foam panels were glued to a non combustible board called "Promatek 
H", density 870 kg/m, with a water based contact adhesive called "Casco 


3880". The non combustible boards were nailed to the light weight concrete 
walls and ceiling before the PUR foam panels were glued. 


Observations during test. 

Time, [muin:s]} Observations 

0:00 Ignition of the burner, 100 kW. 

0:14 The aluminium surface started to get damaged in the burner 
corner, see photo no. 2. 

0:25 Smoke production started. Large portions of the ceiling were 
ignited, see photo no. 3. 

0:38 Flames out the doorway. High SPR was seen. Flame spread 
down the walls, see photo no. 4. 

0:40 _ HRR reached 2 MW. 

1:00 Gas shut off, extinguishment. 

After test: The entire ceiling was charred but not totally consumed. 


Large portions of the aluminium paper surface were still 
intact on the walls, see photo no 6. 


Measured data 

Thickness, mm 41 

Area weight, kg/m? 2:03 

Density, kg/m? 38 (PUR) 

Conditioning ; 


Temperature 20+5 °C 
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Test results, graphs 
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Figure 1 PUR foam panel with al faced paper finish, heat release rate (including the 
ARR from the ignition source), smoke production rate and CO production 
Tate. 
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Figure 2 PUR foam panel with al faced paper finish, ceiling temperatures, wall 


temperatures and heat flux, including the contribution from the ignition 
source. 
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t Photo no 1 Prior to test "PUR foam panel with al faced paper" 


The PUR foam panels were glued to a non combustible board. The non combustible 
boards were nailed to the light weight concrete walls and ceiling. The joints between 
the panels were sealed with aluminium tape. 
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Photo no 2 Time 0:14 "PUR foam panel with al faced paper" 


The aluminium paper surface started to get damaged in the burner corner. 
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Photo no 3 ie 0 cen "PUR foam panel with al faced paper" 


Smoke production started. Large portions of the ceiling were ignited. 


Photo no 4 Time 0:38 "PUR foam panel with al faced paper" 


Flames out the doorway. High HRR and SPR was seen. Flame spread downward the 
walls. 
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Photo no 5 Time 0:48 “PUR foam panel with al faced paper" 


Flash over was reached. Rapid increase in HRR and flame spread. High SPR. 


Photo no 6 After test “PUR foam panel with al faced paper" 


The entire ceiling was charred but not totally consumed. Large portions of the 
aluminium paper surface were still intact on the walls. 
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Test results, ISO 9705 (NT FIRE 025) 


Product 


Mass timber, varnished 


Mounting 


The mass timber was nailed to the lightweight concrete walls and ceiling. 


Observations during test. 


Time, [min:s] 


Observations 


0:00 Ignition of the burner, 100 kW. 

0:25 The lacquer finish had ignited, see photo no. 2. 

0:45 The ceiling tiles above the bummer had ignited. Low SPR. 

1:00 The ceiling was no longer visible due to smoke production, 
see photo no. 4. 

1:30 Downward flame spread was seen on the walls, see photo 
no.Sand6. - 

1:48 Flames out the doorway, see photo no. 7. 

2310 Half of the height of the walls was engulfed in flames. 

2:20 Gas shut off. 7 

2:35 Extinguishment. 

After test: The entire ceiling and slightly more than 50 % of the wall 
surface were charred. Other parts were discoloured or 
undamaged. 

Comments to the graphs 


After approximately one minute the ceiling thermocouple "c2" was not fully in 
contact with the material anymore. After approximately one minute, the reading 
of the heat flux meter at 2 m is higher than the reading of the heat flux meter at 1 
m, mainly because it captured more flame radiation from the burning timber 


than the other one. 
Measured data 


Thickness, mm 


Area weight, kg/m? 
Moisture content, % 


Conditioning 


10 10 
nA 


Temperature 20+5 °C 
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Test results, graphs 


Heat release rate [kW] 


Smoke production rate {m?/s] 


CO production rate [g/s] 


Time {min} 


Figure 1 Mass timber (pine) varnished, heat release rate (including the HRR from the 
ignition source), smoke production rate and carbon monoxide production rate 
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Wall temp. [°C] 


Time [min] 


flux 1m 
— — - flux 2m 


Heat flux [kW/m’} 
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Figure 2 Mass timber (pine) varnished, ceiling temperatures, wall temperatures 
and heat flux, including the contribution from the ignition source. 
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Photo no 1 Prior to test "Mass timber, varnished" 


The mass timber was nailed to the lightweight concrete walls and ceiling. 


Photo no 2 Time 0:25 "Mass timber, varnished" 


The lacquer finish had ignited. 
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Photo no 3 


Time 0:47 
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"Mass timber, varnished" 
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Flames had reached and ignited the ceiling timber above the burner. Low SPR. 


Photo no 4 


SPR had increased. The ceiling was no longer visible. 


Time 1:05 


"Mass timber, varnished" 
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Photo no 5 Time 1:33 "Mass timber, varnished" 


Downward flame spread was seen on the rear wall. 


Photo no 6 Time 1:38 


"Mass timber, varnished" 


Downward flame spread was seen on the left wall. HRR and SPR increased. 
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"Mass timber, varnished" 


Flames were seen out the doorway. Flash over was reached. 
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"Mass timber, varnished" 


Photo no 8 After test 
The room fire was extinguished at 2:35 [min:s]. The entire ceiling and slightly more 


than 50% of the walls were charred. 
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Test results, ISO 9705 (NT FIRE 025) 
Product 

FR chipboard 

Mounting 
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Enclosure 8 


The FR chipboard was nailed to the lightweight concrete walls and ceiling. 


Observations during test. 

Time, [min:s] Observations 

0:00 Ignition of the burner, 100 kW. 

2:00 A smoke gas layer was formed, low HRR, see photo no. 2. 
5-10:00 Low HRR, no flame spread was seen, see photo no. 3. 
10:00 The burner output was increased to 300 kW. 

11:00 HRR and SPR increased. The height of the smoke gas layer 


was approximately 1 meter, see photo no. 4. 


15:00 HRR and SPR had reached maximum levels. Flames were 
seen in the inner half of the ceiling. Limited flame spread on 


the walls, see photo no. 5 


20:00 Gas shut off. 


After test: Afterglow was seen in the material at the burner comer. The 
glowing combustion was extinguished 17 minutes after gas 
shut off. The ceiling was partly charred. Most parts of the 


walls were only discoloured. 
Measured data 
Thickness, mm , 22 
Density, kg/m? 805 
Moisture content, % 6.8 
Conditioning 


Temperature 20+5 °C 


or ott 
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Test results, graphs 


Heat release rate [kW] 


Bes © Sos0s.5 see 


Smoke production rate [m?/s] 


Time [min] 


CO production rate [g/s] 


0 5 10 eS 20 2S 
Time [min] 


Figure 1 FR chipboard, heat release rate (including the HRR from the ignition source), 
smoke production rate and carbon monoxide production rate 
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Photo no 1 Prior to test "FR chipboard" 


The FR chipboard were nailed to the lightweight concrete walls and ceiling. 
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Photo no 2 Time 2:05 "FR chipboard" 


A smoke gas layer was formed, low HRR. 
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Photo no 3 Time 9:57 "FR chipboard" 


Low HRR, no flame spread was seen. Some flames were seen in the material inthe - 
vicinity of the burner. 


Photo no 4 Time 10:47 "FR chipboard" 


The burner heat output had been increased to 300 kW. HRR and SPR increased. The 
height of the smoke gas layer was approximately 1 meter below the ceiling. 
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Photo no 5 Time 14:57 "FR chipboard" 


HRR and SPR had reached maximum levels. Flames were seen in the inner half of 
the ceiling. 


Photo no 6 Time 18:47 "FR chipboard" 


HRR and SPR were still limited. 
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Test results, ISO 9705 (NT FIRE 025) 
Product 

3-layered FR polycarbonate panel 

Mounting 


The polycarbonate panels were screwed to a frame of light steel profiles which 
gave a spacing of 40 mm to the light weight concrete walls and ceiling. 


Observations during test. 


Time, [min:s] Observations 

0:00 Ignition of the burner, 100 kW. 

0:15 The ceiling panel above the burner was deformed. 

1:00 Melted material formed droplets. 

3:00 The polycarbonate burned only in the vicinity of the burner 
flame, see photo no. 3. 

3:40 SPR increased, the ceiling was no longer visible. 

5:00 HRR and SPR decreased, the ceiling was again visible. 

8:00 Some melted material burned on the floor close to the 
burner. | 

10:00 The burner output was increased to 300 kW. Melted material 
moved away from the bumer flame, see photo no 4. 

16:00 Flames were only seen in the burner corner and on the floor 
nearby the burner, see photo no. 5. 

20:00 Gas shut off. 

After test: Some flaming was seen on the floor near the burner. All 


material in the ceiling was melted and had fallen down. 
Large portions of the walls were still intact, see photo no. 6. 


Measured data 

Thickness, mm 16 
Area weight, kg/m? 2.9 
Conditioning 


Temperature 20+5°C 


SP 511 
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Test results, graphs 


Heat release rate [kW] 


0 5 10 15 20 25 
Time [min] 


Smoke production rate [m?/s] 


0 5 10 15 20 Ps) 
Time [min] 


CO production rate [g/s] 


6) S 10 1) 20 25 
Time [min] 


Figure ] 3-layered FR polycarbonate panel, heat release rate (including the HRR from 
the ignition source), smoke production rate and CO production rate 
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Figure 2 3-layered FR polycarbonate panel, ceiling temperatures, wall tempera- 
tures and heat flux, including the contribution from the ignition source. 
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Photo no 1 Prior to test "3-layered FR polycarbonate panel" 


The polycarbonate panels were screwed to a frame of light steel profiles which 
provided for a spacing of 40 mm to the light weight concrete walls and ceiling. 


Photo no 2 Pimest: 30 "3-layered FR polycarbonate panel" 


Melted material formed droplets. Very limited HRR and SPR. 
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Photo no 3 ae 27 4 "3-layered FR polycarbonate panel" 


The polycarbonate panels burnt only in the vicinity of the burner flame. 


Photo no 4 Time 10:20 


"3-layered FR polycarbonate panel" 


The burner heat output had been increased to 300 kW. Melted material mowed away 
from the burner flame. HRR and SPR were still limited. 
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Photo no 5 Time 16:09 "3-layered FR polycarbonate panel" 


Flames were only seen in the burner corner and on the floor near the burner. 
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Photo no 6 After test "3-layered FR polycarbonate panel" 


All material in the ceiling had melted and fallen down. Large portions of the walls 
were still intact. 
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Test results, ISO 9705 (NT FIRE 025) 
Product 

FR expanded polystyrene board, 40 mm 
Mounting 


The polystyrene boards were glued to a non combustible board called “Promatek 
H", density 870 kg/m3, with a water based contact adhesive called "Casco 
3880". The non combustible boards were nailed to the light weight concrete 
walls and ceiling before the polystyrene boards were glued. 


Observations during test. 

Time, [min:s] Observations 

0:00 Ignition of the burner, 100 kW. 

0:20 Burning droplets were formed. 

0:40 Melted material run downwards the walls at the burner 


comer, see photo no. 2. 


1:20 HRR and SPR increased rapidly. The ceiling was no longer 
visible. Melted material was dripping on the floor from the 
entire ceiling, see photo no. 3. Downward flame spread was 


seen. 
teks A few flames were seen out the doorway. 
2:30 The intensity of the fire was decreased. Only flaming in the 


burner corner was seen, see photo no. 4, and a few flames on 
the floor. The ceiling material was all melted. 


10:00 The burner output was increased to 300 kW. 

10:30 HRR and SPR increased. 

11:00 Flaming was seen only in the burner comer, see photo no. 5. 
11:50 A non-combustible backing board fell down from the 


ceiling. Some of the melted polystyrene on the floor was 
then screened from the heat radiation. 


20:00 Gas shut off. 


After test: A few small flames were seen in the material after gas shut 
off. Most material was burned or melted. Some undamaged 
polystyrene slabs were still attached to the lower part of the 


walls. 
Measured data 
Thickness, mm 40 
Density, kg/m? 30 
Conditioning 


Temperature 20+5 °C 
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Test results, graphs 


Heat release rate [kW] 


Time [min] 


Smoke production rate [m?/s) 


CO production rate [g/s] 


Figure 1 FR expanded polystyrene board (40 mm), heat release rate (including the HRR 
from the ignition source), smoke production rate and CO production rate 
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Figure 2 FR expanded polystyrene board (40 mm), ceiling temperatures, wall 


temperatures and heat flux, including the contribution from the 
ignition source. 
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Photo no 3 VE Veale Epa "FR expanded polystyrene board, 40 mm" 


HRR and SPR increased rapidly. The ceiling was no longer visible. Melted material  - 
was dripping to the floor from the entire ceiling. Downward flame spread was seen. 
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Photo no 4 Time 2:31 "FR expanded polystyrene board, 40 mm" 


The intensity of the fire had decreased. Only flaming in the burner corner and a few 
flames on the floor was seen. The ceiling material was all melted. 
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Test results, ISO 9705 (NT FIRE 025) 


Product 
FR expanded polystyrene board, 80 mm 
Mounting 


The polystyrene boards were glued to a non combustible board called "Promatek 
H", density 870 kg/m}, with a water based contact adhesive called "Casco 
3880". The non combustible boards were nailed to the light weight concrete 
walls and ceiling before the polystyrene boards were glued. 


Observations during test. 

Time, [min:s] Observations 

0:00 Ignition of the burner, 100 kW. 

0:15 The material in the ceiling above the burner started to melt, 
see photo no. 2. 

0:30 Melted material run downwards the walls at the burner 
comer, see photo no. 3. 

i25 HRR and SPR increased rapidly. The ceiling was no longer 
visible, see photo no. 4. 

1:45 Melted material was dripping to the floor from the entire 


ceiling. Downward flame spread was seen. A few flames 
were seen out the doorway, see photo no. 5. 


215 The intensity of the fire decreased. 

2:45 Only flaming in the burner comer was seen, see photo no. 6. 

3:00 All ceiling and approximately 50 % of the walls were 
consumed or melted. 

7:00 A small pool fire was seen on the floor nearby the burner 
comer. 

10:00 The burner output was increased to 300 kW. 

10:20 Some increase of HRR and SPR was seen due to burning 

melted material near the burner, see photo no. 7. 

13:00 HRR and SPR started to increase rapidly, see photo no. 8. 

13:21 Flames out the doorway. Downward flame spread was seen 
on the walls, see photo no. 9. 

13:29 Gas shut off, extinguishment. 

After test: Most material had burned or melted, see photo no. 10. 

Measured data 

Thickness, mm 80 

Density, kg/m? 17 

Conditioning 


Temperature 20+5 °C 
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Test results, graphs 


3 
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Heat release rate [kW] 


Smoke production rate [m?/s] 


0 2 4 6 8 10 ib 14 


CO production rate [g/s) 
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Figure 1 FR expanded polystyrene board (80mm), heat release rate (including the HRR 
from the ignition source), smoke production rate and CO production rate 
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Photo no 1 Prior to test "FR expanded polystyrene board, 80 mm" 


The polystyrene boards were glued to a non combustible board. The non combustible 
boards were nailed to the light weight concrete walls and ceiling. 
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Photo no 2 Time 0:16 "FR expanded polystyrene board, 80 mm" 


The material in the ceiling above the burner started to melt. 


Yeh REP OMMESaE aOR: Cite ae 


te fond 


WAbertors-----—- 6-4 


ig oapacebhagd oust et 
tal 2) of Baim! Siw a 
J 


es 


a pid Wee » 


¥ 
r 
-* 


as an fr 
2 ; 


mae goat 9 


oe ioe Sa 


< 
7) 


"rom O8 .bisod ecanonitot bebaagts JEP". ws . w: f T 
Figure FE cxpanee s polysegre Linon), Bratt 
‘re 


Mi ri 
roan av “ener . ) 


DatunV Date Beteckning/Reference Sida/Page 


1996-01-26 95R22049 5 (8) 
Enclosure 11 


Photo no 3 Time 0:36 "FR expanded polystyrene board, 80 mm" 


Melted material run downwards the walls at the burner corner. 
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Photo no 4 Time 1:26 "FR expanded polystyrene board, 80 mm" 


HRR and SPR increased rapidly. The ceiling was no longer visible. 
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Photo no 5 Varied 5.1 "FR expanded polystyrene board, 80 mm" 


Melted material was dripping to the floor from the entire ceiling. Downward flame 
spread was seen. A few flames were seen out the doorway. 
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Photo no 6 Time 2:46 "FR expanded polystyrene board, 80 mm" 


Only flaming in the burner corner was seen. 
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Photo no 7 


Time 10:20 


"FR expanded polystyrene board, 80 mm" 


The burner heat output had been increased to 300 kW. Some increase of HRR and 
SPR was seen due to melted material near the burner that ignited 


Photo no 8 


Time 13:11 “FR expanded polystyrene board, 80 mm" 


HRR and SPR started to increase rapidly. 
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Photo no 9 Pbumes!3; 23 "FR expanded polystyrene board, 80 mm" 


Flash over was reached. High HRR and SPR, flames out the doorway. Downward 
flame spread was seen on the walls 
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Photo no 10 After test "FR expanded polystyrene board, 80 mm" 


Most material had burnt or melted. 
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Test results, ISO 9705 (NT FIRE 025) 

Product 

Plywood 

Mounting 

The plywood was nailed to the lightweight concrete walls and ceiling. 

Observations during test. 

Time, [min:s] Observations 

0:00 Ignition of the burner, 100 kW. 

0:45 The material in the burner comer and in the ceiling above 
the burner was ignited, see photo no. 2. 

1:30 More than 50 % of the ceiling was ignited and flame spread 
was starting downward the walls, see photo no. 3. 

2:14 Flames were seen out the doorway. SPR increased. Flame 
spread was seen approximately 1 m down the walls, see 
photo no. 4. 

2:58 The ceiling and most parts of the walls were engulfed in 
flames. HRR peaked at approximately 2 MW. 

3:00 Gas shut off. 

3:05 Extinguishment. 

After test: Most parts of the wails and ceiling were charred. 

Measured data 

Thickness, mm 14 

Density, kg/m? 440 

Moisture content, % 11.3 

Conditioning 


Temperature 20+5°C 
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Test results, graphs 
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Heat release rate (kW) 
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CO production rate [g/s] 
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Figure I Plywood, heat release rate (including the HRR from the ignition source), 
smoke production rate and carbon monoxide production rate 
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Wall temp. [°C] 
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Figure 2 Plywood, ceiling temperatures, wall temperatures and heat flux, 
including the contribution from the ignition source. 
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Photo no 1 Prior to test "Plywood" 


The plywood was nailed to the light weight concrete walls and ceiling. . 


Photo no 2 Time 0:46 


"Plywood" 


The material in the burner corner and in the ceiling above the burner had ignited 
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Photo no 3 Aime .25 "Plywood" 


More than 50 % of the ceiling had ignited and flame spread was starting downward 


the walls. 


oes 


Photo no 4 
Flames were seen out the doorway. SPR increased. Flame spread was seen 
approximately 1 m down the walls. 


Time 2:14 "Plywood" 
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Photo no 5 ime 2:58 "Plywood" 


Flash over was reached. The ceiling and most parts of the walls 
flames. HRR peaked at approximately 2 MW. 


Photo no 6 After test : 


Most parts of the walls and ceiling were charred. 
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Test results, ISO 9705 (NT FIRE 025) 
Product 

FR Plywood 

Mounting 


The FR plywood was nailed to the lightweight concrete walls and ceiling. 


Observations during test. 

Time, [min:s] Observations 

0:00 Ignition of the burner, 100 kW. 

0-10:00 Only limited buming in the burner corner was seen. Low 
SPR. 

10:00 The burner output was increased to 300 kW. 

10:10 Flame spread in the ceiling was seen. SPR increased, see 
photo no. 3. 

10:35 Flames out the doorway. 

10:45 Downward flame spread was seen on the walls. SPR 
decreased, see photo no. 4. 

11:10 Gas shut off. 

11:16 HRR decreased, see photo no. 5. 

£1:53 The material started to self extinguish, see photo no. 6. 

After test: Most parts of the walls and ceiling were charred. 

Measured data 

Thickness, mm 15 

Density, kg/m? 460 

Moisture content, % 9.8 

Conditioning 


Temperature 20+5°C 
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Test results, graphs 


Heat release rate [kW] 


W 
8 


—_— —_ 
o NS 


Smoke production rate [m?/s] 


CO production rate [g/s) 


0 Pe 4 6 8 10 fe 
Time [min] 


Figure 1 FR plywood, heat release rate (including the HRR from the ignition source), 
smoke production rate and carbon monoxide production rate 


11 


Sf 


PR RIS DENRADT 
’ ¢ 25> 29 ss se = DatunvDate Beteckning/Reference Sida/Page 

Loos oO 1996-01-26 95 R22049 4 (6) 
° _ Enclosure 13 
< Zz 
es x 

° & 

VG-Fro* 


Photo no 1 Prior to test "FR Plywood" 


The FR plywood was nailed to the light weight concrete walls and ceiling. 
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Photo no 2 hime [313 "FR Plywood" 


Limited burning in the vicinity of the burner flame. 
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Photo no 3 Time 10:12 "FR Plywood" 


The burner heat output had been increased to 300 kW. Flame spread in the ceiling 
was seen. SPR increased. 
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Photo no 4 Time 10:45 "FR Plywood" 


Flash over was reached. Downward flame spread was seen on the walls. SPR 
decreased. 
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Photo no 5 Time 11:16 


"FR Plywood" 
The gas burner had been shut off due to flash over. HRR decreased. 
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Photo no 6 "FR Plywood" 


The material self-extinguished. ~ 
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Appendix D — Quintiere’s Fire Growth Model 


D.1 Quintiere, J. G., “A Simulation Model for 
Fire Growth on Materials Subject to a 
Room/Corner Test”, Fire Safety Journal, 
Volume 18, 1992. 

D.2 Fire Growth Model Source Code 
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D.1 — Quintiere J G., “A Simulation Model for Fire Growth on 
Materials Subject to a Room/Corner Test”, Fire Safety 
Journal, Volume 18, 1992. 
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ABSTRACT 


A mathematical model has been developed to simulate fire growth on 
wall and ceiling materials when subject to a room-corner fire test 
exposure. The model predicts the area of burning, the upper layer gas 
temperature, and the rate of energy release as a function of time. 
Material fire property data are developed from apparatuses described in 
ASTM E 1321 and E 1354. The results compare favorably to experi- 
mental data generated in Sweden for 13 materials tested. Furthermore, 
the model shows the sensitivity to flashover’ for thin materials relative 
to small variations in their property data. 


NOTATION 


area 
parameter defined in eqn (14) 
specific heat 

depth of room 

side of square burner 

acceleration due to gravity 
convective heat transfer coefficient 
height of room, vent 

thermal conductivity 

empirical constant, eqn (7) 
effective heat of gasification 
empirical power, eqn (7) 
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q heat 

Q energy release 

t time 

fe temperature 

w width of room 

x lateral position 

y upward position z 
z downward position ° 
p density ; 
T dummy variable for time, eqn (1) “ 
AH heat of combustion 

Subscripts 

b burn-out 

if flame 

ig ignitor, ignition 

min minimum 

p pyrolysis 

S surface 

S, 0 surface responding to ignitor flame heat flux 

O initial 

co ambient 

Superscripts 

(‘) per unit time 

( )’ per unit width 

(FF per unit area 


INTRODUCTION 


For many years the regulation of interior finish materials applied to the 
walls, floor and ceiling of a building has been determined by a 
flammability test. The criterion for acceptance is usually based on a 
rating scale for the flammability test. This gives a relative level of & 
performance for the material’s potential for fire growth. Other fire — 
safety elements such as smoke obscuration and toxicity are treated in a = 
similar fashion. The limitations of a relative ranking scale based on a s 
single test is that it is not necessarily applicable to the material’s fire 
growth potential in its actual application. 

Lack of confidence in the test ranking approach led to the develop- 
ment of a standard room test to investigate the material’s fire growth 
potential." Although the room test is conducted at a more realistic 
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scale, it still only offers one result, e.g. the time to achieve a specific 
energy release or flashover. The effect of the ignition source and room 
geometry are not known without conducting a series of expensive tests. 

Consequently, an alternative to conducting full-scale tests has the 
hope of using material fire property data to predict such results for a 
wide range of conditions. The first attempt for making predictions of 
room fire growth over interior finish materials was done by Smith.* He 
employed material energy release rate data from a dynamic calorimetry 
apparatus.* More recently Karlsson and Magnusson‘ developed a model 
for room fire growth on wall and ceiling materials which incorporated 
data from two new fire property test methods.°*° Both of these models 
used empirical methods for computing flame spread. As an attempt to 
improve this modeling approach, Cleary and Quintiere’ presented a 
simple, but complete accounting for all modes of flame spread which 
govern growth on a wall and a ceiling. Although successful agreement 
was found with data, their model lacked a direct accounting of room 
thermal feedback, and selected energy release rate data from the Cone 
Calorimeter at an arbitrary irradiance level. The simulation model 
described in this paper attempts to eliminate these two limitations. 
However, as will be seen, the new model requires the specification of 
the ignition source and spreading flame heat flux which will still 
introduce some uncertainty. Although Smith? computed the flame heat 
flux in his model, this author does not believe that accurate results are 
possible by purely theoretical means. Instead, heat flux correlations 
developed from controlled experiments will be more reliable. 

The validity of any model must be judged on the soundness of its 
components, its completeness, and its ability to predict experimental 
results. Two of the models described above*’ have been evaluated 
against an extensively instrumented series of room-corner fire experi- 
ments for 13 materials.* A complete set of fire property data was also 
available for the 13 materials.’? Both of the models,*’ as well as 
empirical models by Wickstrém and Goransson’® and Ostman and 
Nussbaum,"' have all given good predictive results for the time to 
flashover in these experiments. All of these models have used the rate 
of energy release per unit area (Q”) as a significant input parameter. 
The issues that remain are (1) what is the most appropriate form of the 
fire property data, (2) what is the complete set of properties required, 
and (3) what predictive model is the most general. 

An extension of the model presented by Cleary and Quintiere’ will 
be described in this paper. It will demonstrate the way fire property 
data are required from ASTME 1321 (ISO WD 5658 Part II) and 
E 1354 (ISO IS 5660), and will include the effect of thermal feedback 
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due to the increase of temperature within the room. These results will 
also be compared to the Swedish room-corner fire experiments.® 


DESCRIPTION OF THE MODEL 


The model simulates the ignition, flame spread, burn-out, and burning 
rate of wall and ceiling materials subject to a corner fire ignition source 
in a room. The flame pyrolysis and burn-out fronts are computed with 
respect to two modes of flame spread. One mode includes upward 
spread, spread along the ceiling, and spread along the wall-ceiling jet 
region. At this time, no distinction for these different configurations is 
made in the model and they are universally treated as governed by 
upward flame spread. The second mode of spread is composed of 
lateral spread along the wall and subsequent downward spread from the 
ceiling jet. Again, the same relationship will be considered for both. In 
this fashion, the pyrolysis and burn-out areas are computed. 

The energy release rate per unit area is considered as a function of 
time. The energy release rate per unit area is governed by the flame 
heat flux and the radiative feedback from the heated room. Flame heat 
flux is considered uniform over the pyrolysis area, and uniform over the 
extended flame length. Two values are selected: 60 kW/m? over the 
pyrolysis area and associated with the square burner corner ignition 
flame, and 30kW/m7? for the extended flame spread which governs 
upward flame spread. 

The room thermal feedback controls both the rate of spread through 
a computation of the material surface temperature ahead of the flame, 
and the rate of energy release per unit area through radiative heat 
transfer from the gas layer in the room. Global models are considered 
for room surface temperature and gas layer temperature. The radiative 
effects are considered to be maximized to give an upper limit for its 
effect. It appears that the thermal feedback effect is not significant 
compared to the flame heating effects until conditions representative of 
the onset of flashover, e.g. a gas temperature of 500°C and a 
corresponding blackbody irradiance of 20kW/m*. Hence a more 
detailed representation of the room thermal feedback may not be 
productive at this time. 

The details for each component of the analysis will be described in 
the following sections. 


FIRE SCENARIO 


The corner fire scenario considered in this analysis is based on tests 
described by Sundstrém.* A material lined the walls and ceiling of a 
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room 2-4 m X 3-6 m X 2-4m high with a doorway opening 0-8 m x 2-0 m 
high. A square burner, 0-17m on a side, located at the floor in a 
comer, subjected the wall to an ignition flame. The test scenario 
prescribed a 100 kW ignition flame for 10 min, followed by a 300 kW 
ignition flame. Each of 13 materials was examined up to flashover or 
burn-out under this ignition scenario. Flashover was experimentally 
identified as coinciding with 1 MW energy release rate from the room, 
and this will correspond to the flashover time computed in the situation. 


IGNITION BURNER 


The heights of the burner flame corresponding to 100 kW and 300 kW 
were modified from the process analysis by Cleary and Quintiere’ to be 
1-3 and 3-6 m respectively. The latter value corresponds to a correlation 
by Hasemi and Tokunaga” for corner burner fires. Their correlation 
gives a flame tip height of 1-9 m for the 100 kW fire compared with the 
1-3m value retained in the current simulation. A later analysis by 
Hasemi and Tokunaga’ gives 2‘8m and 5-9m for the flame tips and 
2:1m and 4-4m for the continuous flame height corresponding to 
100kW and 300kW for the 0-17m square corner burner. These 
variations indicate the uncertainty, and suggest the need for additional 
study. Moreover, in the current situation, the flame length will be 
assumed to be representative of a vertical wall flame. The influence of a 
ceiling results in a flame extension in the model equal to the distance 
the representative wall flame is above the physical ceiling. 

The burner is assumed to initially prescribe a uniform heat flux to the 
wall over a region defined by the flame height and the width of the 
burner. Based on a study by Williamson et al.,'* the burner heat flux 
has been taken as 60 kW/m7’. In general, it appears that the burner 
flame heat flux depends on the size of the burner, the flame height or 
energy release rate, and the type of fuel supplied. Hence, a general 
simulation of room corner fires must be able to represent these effects. 
Moreover, the heat flux from the flames due to the burning surfaces in — 
the corner and along the ceiling must also be predictable. At this time, 
only rational assumptions can be made for these heat fluxes. 


IGNITION SIMULATION 


The simulation model assumes a burner heat flux prescribed over a 
region of width (D), the burner side, and a flame height of 1-3 m for 
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the 100kW ignitor. The surface temperature of the material is 
computed by” 


1’ q(t) 
Pi, gle | dt 1 
Vakpc Jo Vt—T (1) 


where q(t)=4i,+ 0(T*— T?.), Gi, is the ignitor flame heat flux 
assumed at 60 kW/m? and T is the temperature of the upper gas layer in 
the room. 

Blackbody radiation has been assumed with a configuration factor 
equal to one to maximize the room feedback effect. When T,., equals 
the ignition temperature (7;,), the time for ignition by the burner (¢,,) 
is found. At this time, the material begins to contribute energy and the 
flame spread process commences. Figure 1 indicates the path of flame 
extension from the burner and burning material. The ceiling jet region 
is approximated at a depth of 0-08H where H is the ceiling height. This 


Quadrant of ceiling 


d= 3.6m 


w=2.4m 


Doorway 


2.0 m 


0.8 m 
D=0.17m 


Room, burner and doorway dimensions 


Fig. 1. Room and burner configuration. Dashed lines indicate path of upward flame 
extension. 
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is derived from Alpert’s study of axisymmetric ceiling jets.’® It is 
approximated as a constant depth because it varies slowly with distance 
from the fire source. 


UPPER LAYER GAS TEMPERATURE 


In general, the room gas and surface temperature vary with position 
and time, and are highly coupled with the energy release rate of the 
growing fire. Since it is expected that room thermal effects are not 
significant during early fire growth, the simplest representation for 
upper layer gas temperature is used. This eliminates the need for a 
comprehensive analysis by compartment zone or field models. How- 
ever, the other aspects of the simulation model for fire growth could 
ultimately be incorporated into more comprehensive compartment fire 
models. The approach taken here is based on the compartment 
temperature correlation of McCaffrey et al.'’ with the coefficient 
modified for corner fires as done by Karlsson and Magnusson.* The 
correlation is given as 


QO kpc/t A, ie (2) 


2/3 
fe ae mg ce 
p.CVB TAG Hi, p.CpVE A, A, 


where 


QO is the total energy release rate 

A, is the room surface area | 

A, is the area of the opening 

H, is the height of the opening 

kpc is the thermal inertia of the room lining materials 

pc,Veg is 3-44 kW/m>?— 

C is the coefficient taken as 2-2 for these corner fires (compared to 
1-63 for room-centered fires) 


The higher value of C is expected for corner fires compared to centered 
fires because of the lower entrainment rate of air into the corner fire. In 
other words, corner fires are hotter. 


ENERGY RELEASE RATE 


The total energy release rate is composed of the energy from the 
ignition burner and the energy from the wall and ceiling materials. This 
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is given by ) 
O(t) ia On os O"A,(t) 


where Q,, is the ignition burner energy release rate, O”(t) is the energy 
release per unit area of the material and A, is the pyrolysis area. 


ENERGY FROM MATERIAL 


The rate of energy release rate per unit area is assumed constant at any 
instant of time and uniform over the pyrolysis area. It is dependent on 
the net heat flux incident on the material which is composed of the 
flame heat flux, room thermal feedback flux, and the re-radiation loss 
assumed to occur at the ignition temperature of the material. The 
relationship for Q” is developed in terms of peak values of Q” found 
from the Cone Calorimeter for different irradiance levels (incident 
external radiant heat flux) as shown in Fig. 2. Peak values for Q” are 
denoted on the figure with some suggested degree of uncertainty. In 
~~Fig. 3, these peak values are plotted against the Cone irradiance level. 
The linear plot suggests that the net flame heat flux in the Cone 
Calorimeter is a constant. From the theory of burning rate, the slope is 
AH/L where AH is the heat of combustion and L is taken to be an 
effective heat of gasification. For the steady burning rate of liquid fuels, 


800 Textile wall covering on gypsum board 
Irradiance Q”"=48x103 —8:3x103 kJ/m2 

rn 75 kW/ m2 
600 | 

' 

rh 

_ 2 

400 ae 50 kW/m 


200 


Energy release rate per unit area (kW/m?) 


oO 100 200 300 400 
Time (s) 


Fig. 2. Example of Cone Calorimeter data taken from Ref. 9. 
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700 


Textile wall -covering 
on gypsum board 


600 


“Siope, AH/L 
AH =13 kJ/g 
L=15-19kJ/g 


(kW/m?) 
A 
re) 
(eo) 


Peak energy release rate per unit area 


fe) 20 40 60 80 100 
External heat flux (kW/m?) 


Fig. 3. Example of linear relationship between peak OQ” and irradiance for Cone 
Calorimeter data. : 


L is a thermodynamic material property. Here, L is a modelling 
parameter that enables the determination of peak values of Q”. Since 
the heat of combustion can be derived independently from the Cone 
Calorimeter by dividing the instantaneous energy release rate by the 
mass loss rate, the slope of the data in Fig. 3 can be used to determine 
L. In that example for the textile wall-covering on gypsum board, L is 
found to be 1-5 kJ/g for the absolute peak values of Q”, but can range 
as high as 1-9 kJ/g if lower peak values are used. 

If it is assumed that AH/L is an effective material property 
characteristic of the burning behavior which enables the computation of 
the peak values under all heat flux conditions, then 


©, ANN 
Q"= WE (q;— oT;, + oT“) (4) 


where gq; is the incident flame heat flux over the pyrolysis region (a 
specied constant), oT;, is the re-radiation flux loss, and oT% is the 
incident heat flux from the room, maximized as a blackbody gas with a 
configuration factor of one for all locations. As the compartment 
temperature, 7, varies with time, the value for QO” varies with time. 
Since the material initially burns due to the ignition burner flame, 
and the material flame heat flux is likely to be similar to the burner in 
its vicinity, 60kW/m? is also assumed for gq; here. As propagation 
moves far beyond the burner flame region, or to the ceiling or ceiling 
jet regions, some variations in q; are likely. Also this may depend on 
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the extent of the pyrolysis region. Such information on flame heat flux 
is not currently available. Hence, burning is assumed to occur based on 
the ignition burner flame heat flux for the entire pyrolysis region. 


THE PYROLYSIS AREA, A, 


The pyrolysis area is computed from the configuration of the pyrolysis 
and burn-out fronts. Figure 4 illustrates a construction of the pyrolysis 
area as bounded between the pyrolysis and burn-out fronts. The 
pyrolysis front represents the position where pyrolysis has just begun; 
and the burn-out front represents the position where burn-out has just 
occurred. The initial pyrolysis area occurs at t=4,, for the region 
bounded by y,,, = 1-3 m, the burner flame height, and x,,, = 0-17 m, the 


pedcave 


fa 


Fig. 4. Illustration of the pyrolysis and burn-out regions constructed from the 
pyrolysis and burn-out fronts. 
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dimension of the burner side. For the corner, this area is multiplied by 
2 to account for the two symmetric walls. As spread and burn-out 
continues, Fig. 4 illustrates the upward spread pyrolysis front, y,, to be 
greater than the room height. Consequently, it has been extended along 
the ceiling by a radius y, — H, and along the wall—ceiling jet region by 
the same distance. The burn-out front associated with the upward 
spread component, y,, is still less than the room height for this case. 
Corresponding pyrolysis and burn-out fronts for lateral spread are 
denoted by x, and x, respectively in Fig. 4. A downward pyrolysis 
front, z, is also shown. It commences when y, = H. The pyrolysis and 
burnout areas are formed by connecting the y-fronts to the x- and 
z-fronts by straight lines as shown in Fig. 4. The pyrolysis area is 
computed for the case shown in Fig. 4 as follows: 


Ap= Ap, — Ap. + Apy, + Apc; (5) 
where 


Apia cliixoertatx, — <p Veo + (4)(H —y. \(Xp= Xp0)} 
Apz = 2[YpXpo + (Xe — Xp.0)Yp.0 + (2) (Yo — Yp.o) (Xb — Xp.o)] 
Apn = 2[(% — 1)(0-08H) + (2)(2p)(¥» — H) 

— (3)(0-08H + 2,)’(Xp — Xp.0)/(H — yp.0)] 


where the last term of Ap;, is an approximation to the tnangular 
overlap region 


Apc, = min [(77/4)(y, — H)’, wd] 


where wd is the area of the ceiling. Note that Ap, represents the 
burn-out area, Ap;, is the wall pyrolysis area associated with the ceiling 
jet and Apc, is the pyrolysis area for the ceiling. 

The pyrolysis region can have various configurations depending on 
the values of the fronts. The complete set of configurations is given in 
the Appendix. 

The differential equations governing the fronts will now be examined. 


UPWARD FLAME SPREAD 


The pyrolysis front y, is applied to upward spread along the wall, the 
wall—ceiling jet, and the ceiling. It is assumed that a universal theory of 
upward flame spread governs all of these regions. As more sophistica- 
tion is developed, more appropriate flame spread relationships could be 
applied. The governing equation is based on a constant flame heat flux 
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applied over the flame extension region beyond y, whose initial 
temperature is 7,."° . 


dy; imp 
da > hi (6) 
where 
X T,. — T.7° 
t =—k 8 | 
gs 4 pe gi 
and 
k y" + 9 wif n < yin 
Ys = Yo +{ ey or at ig, Soe (7) 
k{Q (= Wa)ie YW =k Qiz 


T, is the global room average surface temperature computed by eqn (1), 
but with 


g(t) Ol een a! (8) 


and h, = 0-01 kW/m’°K as the convective heat transfer coefficient. 

g; is the flame heat flux beyond the burning region. For vertical wall 
flame spread, this heat flux has been found to be 25+5kW/m’ for 
flame extensions less than 2m, and is relatively independent of the 
fuel.’® Based on this, g; was specified as 30 kW/m? in the simulation. 

y, is the flame length in the upward or wind-aided direction. Qj, is the 
energy release rate for the burner which is equivalent to a line-source. 
It is determined, based on flame length, such that the burner flame 
length corresponding to Q,, is equal to k,Q;". The two possibilities for 
y, indicated above either include the burner effect if y, is less than or 
equal to the burner flame length, or it does not include the burner 
effect. In the latter case, the burner flame and the flame spreading on 
the wall material are not contiguous. The flame length for wall flames is 
given such that k,=0-067 (m°/kW’)'? and n=, or approximately 
k,=0-01 m?/kW and n = 1.'®-!° 


UPWARD BURN-OUT POSITION 
The position of the burn-out front can be approximated by the 


difference equation 


dy, LB Yo(t + te) = Yo(t) 


dt tp oa) 


but it can be shown that y,(t + ¢,) is identical to y,(t) since burn-out will 
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occur at that position in the time interval t,. Therefore 


dy, Y(t) — y(t) 
— = 9b 
dt te isp) 


gives the differential equation for the burn-out front where 
t, = OF ae 
and Q” is the total available energy per unit area which is assumed 


constant for a given material. It is determined by the area under the 
curves of the Cone Calorimeter data shown in Fig. 2. 


LATERAL OR DOWNWARD SPREAD 


The lateral pyrolysis position is determined by 
St aed a 
dt kpc(T, — T,)’ 
where ® and 7, nin are material dependent properties derived from the 
test procedure of ASTM E-1321.° 


The downward pyrolysis position is given for t>t,, the time when 
Yp = H as 


for T, = T; min (10) 


zene) — Xp(tu) (11) 


LATERAL OR DOWNWARD BURN-OUT 


Based on the same logic as in developing eqn (9b), 
Dee Mery Xb 


12 
dt te at?) 


and the downward burn-out front is given, as in eqn (11), by 
Zp = X(t) — Xo(tn), (13) 
where f;, is the time when y, = H. 


INITIAL CONDITIONS FOR THE FRONTS 


For 0S tai: 


yp = Xp 
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At t=t,,. (ignition due to burner flame): 


Yp = Yp.o = 13m 


For 0S ¢t < (tigo + bb,0): 
Y =X, =0 


At t = (tigo + tv,0) (initial burning region extinguishes) 
where 


bo.o = OHO" (tie 0) 


Xp = Xpo 

Yo = Yp.o 
For OStSt, 

zZ>=0 
and for 0=t =f), 

2 =90 


These conditions constitute the values of the fronts over initial time 
periods, and the initial conditions for the differential equations. 


SOLUTION METHODOLOGY 


The equations to be solved constitute four ordinary differential equa- 
tions, one integral equation, and one algebraic equation. They are 
summarized as follows: 


dy 

eqn (6): 4° =Sfi(p, Ye T, T) 
d 

eqn (9b): = fips Yes T) 
dx 

eqn (10): > *=(T.) 


dx 
eqn (12): an = fi(Xp, Xv, T) 


eqn (2): T = f(T, Yp> Yo» Xp» Xb, t) 
eqn (1), (8): T=f,(7,, T, t) 
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First eqns (1) and (2) are solved until the burner ignites the wall 
material. Equation (1), the integral equation is integrated over the time 
steps using the Trapezoidal Rule, and a Gauss-Siedel iterative process 
is employed to obtain the new value of 7;. A Regula Falsi iterative 
method was required for eqn (2) to obtain consistent convergence for 
We 

Once ignition occurs the differential equations are integrated by a 
second order Runga—Kutta method, and the entire set are solved 
simultaneously advancing in time. Estimated values for T and T; are 
used in determining the new values of the x and y variables, then these 
values are used to compute the new values for Tand T.. | 

In the simulation of the tests, the calculation is continued for 1000s 
or until the total energy release rate reaches 2MW. At 600s, the 
burner ignition source is increased from 100 to 300 kW as prescribed in 
the scenario. ; 


MATERIAL PROPERTY DATA 


The material property data needed are displayed in Table 1 for the 13 
materials tested. Most of these data were computed from various 


TABLE 1 
Flame Spread and Heat Release Properties of Swedish Fire Test Materials* 


Material re kpc ? T. AH. ja Q” 


4 S.muin 


(°C) (kW/m*K)* (kW"/m*) (TC) — (Ke) (Kg) (W/m) 
Particle board (PB) 405 0-626 8 180 14 5-4 21-2x 10° 


Insulating fiberboard (IFB) 381 0-229 14 90 14 4-2 26-8x 10° 

Medium density fiberboard 361 0-732 11 80 14 4-2 =10° 
(MDFB) 

Wood panel (spruce) (WPS) 389 0-569 24 155 15 6-3 21-2x 10° 

Melamine covered particle 483 0-804 <1 435 11 48 =6-0x 10° 
board (MELPB) 

Paper covered gypsum board 388 0-593 0-5 300 10 4-8 72x 10° 
(PAPGB) 

PVC covered gypsum board 410 0-208 25 300 13 3-7 4-6 x 10° 
(PVCGB) 

Textile covered gypsum 406 0-570 9 270 13 1-5 8-3x 10° 
board (TEXGB) 

Textile covered mineral 391 0-183 6 174 25 2:8 9-3x 10° 
wool (TEXMW) 

Paper covered particle board 426 0-680 13 250 13 6°5 =10° 
(PAPPB) 

Polyurethane foam (ngid) 393 0-031 3 105 13 3-1 1-4 x 10° 
(PU) 

Expanded polystyrene (EPS) 482 0-464 31° ~130° 28 15 3-2 x 10° 

Gypsum board (GP) 469 0-515 14 380 ” 4-8 2:8 x 10° 


“ From Refs 7 and 9 except as noted. 
° Based on Cone data at 50 kW/ m7” irradiance. 
“ From Ref. 20. 
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sources by Cleary and Quintiere.’ The current analysis developed the 
properties L and Q” by analyzing the Cone Calorimeter data of 
Tsantaridis and Ostman? as discussed relative to Figs 2 and 3. The 
accuracy of L and Q” are limited due to the few available data points and 
the coarse approximation for integrating the curves. Consequently, a 
sensitivity analysis was conducted for some materials by varying L and 
Q” from the ‘base-line’ values in Table 1, e.g. 0-75 L implies a value 
0-75 of the baseline for L. It will be seen that thin materials, which 
have burning times of the order of 10s, or Q” values of 10° to 
10* kJ/m?, are very sensitive to variations in Q” and L with respect to 
the flashover time produced. 


TABLE 2 
Time to Achieve a 1 MW Fire for the Swedish Room Tests 
Material Experi- Base-line 1-25L  1:250/._.J-25L/* "0-756 
mental (s) (s) 0-750” 0:5Q” (s) 
time (s) (5) 
(s) 
Particle board (PB) 157 121 167 88 
Insulating fiberboard oe) 29 36 
(IFB) 
Medium density 131 91 120 
fiberboard (MDFB) 
Wood panel (spruce) 131 110 156 
(WPS) 
Melamine covered 465 Z22 295 
particle board 
. (MELPB) : 
Paper covered gypsum 640 613 620 622 625 
board (PAPGB) 
PVC covered gypsum 611 30 38 42 602 
board (PVCGB) 
Textile covered gypsum 639 41 47 52 606 
board (TEXGB) 
Textile covered mineral 43 12 14 14 14 
wool (TEXMW) 
Paper covered particle 143 222 346 148 
board (PAPPB) 
Rigid polyurethane 6 4 4 
foam (PU) 
Expanded polystyrene 115 Ag 47 
(EPS) 
Gypsum board (GB) » 642 649 650 726 kW? 
650 s 


* Did not reach 1 MW. 
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The property results tabulated in Table 1 were derived from several 
sources and in some cases data from two sources did not agree, hence 
decisions were made in favor of the most consistent data.’ Sufficient 
data were not available for expanded polystyrene (EPS), so generic 
data were taken for @® and T,,,,, from another source for 
completeness.” However, those data were for horizontal flame spread, 
and the unusually high ® value may be due to a melting effect. It is not 
clear that this ® would apply to the vertical case. 


RESULTS 


Simulations were run for all thirteen materials for the room-corner fire 
tests scenario.* Sensitivity analyses for variations in Q” and L were 
performed for some materials. The results are tabulated in Table 2 and 
plotted in Fig. 5. Thick homogeneous materials yield predicted flash- 
over times (times to reach 1 MW) generally lower but within 50% of the 
experimental values. For the variables in L and Q” used (+25%), the 
predicted accuracy is still within 50%. The Melamine covered Particle 
Board (MELPB) burned in two modes due to the laminated construction 
and has a large time difference between the predicted and experimental 
values. The uniform burning assumption used in the model may be too 
crude for the actual burning variations of this combustible composite. 


A Base line a 
way + 125 éh 
© 125L/0:75Q" fa 
600 0 1251/0500" Pt ia 
X O75L Aros 
ea 
500 2 
3 o i 
‘ve 400 a 4 oO 
E x 4 
- + ‘a w Gt 
% 300 Y« re O| |x 
~ “ff >| |W 
Vv : a) |e 
Street wont apo 
£ 200 % 


fe) 100 200 300 400 500 600 fers 
Experimental time (s) 


Fig. 5. A comparison of the predicted flashover times with the experimental times to 
reach 1 MW. 
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The Paper covered Particle Board (PAPPB) is another example of a 
combustible composite that presents a challenge in modeling its burning 
behavior. The expanded polystyrene melts and spills over the floor 
which influences the flashover time, and the melting effect is not 
addressed at all in the model. 

The thin materials on inert substrates display an interesting effect. In 
most cases, flashover did not occur until after the 300 kW burner level 
was initiated in the test. However, the simulated results for PVC 
covered Gypsum Board, Textile covered Gypsum Board, and Gypsum 
Board are only consistent with the experimental results for 1-25L and 
0-5Q” input data as seen in Table 2. Other values of a" and L yielded 
earlier flashover times. 

These variations in flashover times associated with relatively small 
changes in the input data are due to the tendency of upward flame 
spread to accelerate or decelerate. From the linearized theory discussed 
by Cleary and Quintiere,’ the parameter 


b= k,Q" -1- tig! te (14) 


controls acceleration, leading to acceleration if b>0O and ultimate 
extinction if b <Q. In the current simulation, b varies with time and for 
b near zero, small perturbations in the input data can lead to big 
differences in fire growth. A material having a b near zero may not 
appear to give repetitively consistent results even in actual full-scale fire 
testing. 

It should be noted that for these materials studied, fire growth due to 
lateral spread was not a significant factor until the onset of flashover. 
The relatively low flame spread rate, and the need to achieve a 
minimum surface temperature (7, min) before spread could commence 
are reasons for its lack of significance in this simulation model. It may 
be argued that an improved radiative exchange model for the corner 
fire might lead to more significant lateral flame spread effects, but 
enhanced radiation due to the corner flames is countered by the 
maximized gas layer radiation assumed in the model. | 

Figure 6a through 6d show the nature of fire growth for the Textile 
on Gypsum Board for case 0-5Q”/1-25L. In Fig. 6c, the fire appears 
headed for extinction, but when the 300 kW burmer fire is initiated at 
600s, its flame extends beyond the pyrolysis region in the ceiling jet. 
This causes acceleration of the pyrolysis front and flashover. Figure 7 
shows the corresponding predicted rate of energy release compared to 
both the base-line data input and the experimental results. The 
predicted time for ignition and initial rate of spread are faster than the 
experimental results. The sensitivity to flashover with Q” and L is 
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Fig. 7. Comparison of predicted energy release rates with experiments results for 
textile covered gypsum board. 


apparent. This is also seen for Gypsum Board in Fig. 8. The base-line 
case leads to flashover while the 0-5Q”/1-25L case only achieves about 
750 kW compared with 500 kW for the experimental result. 


CONCLUSIONS 
A simulation model has been developed that successfully addresses fire 


growth on walls and ceiling materials in a room corner-fire scenario. All 
modes of fire spread are modeled approximately, but validated ceiling 
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Fig. 8. Comparison of predicted energy release rates with experimental results for 
gypsum board. 
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and wall—ceiling jet flame spread models need to be developed. 
Assumptions have been made based on limited information for flame 
heat transfer rates to cause burning and spread, and more complete 
experimental results are needed. A technique has been included to use 
Cone Calorimeter data in the model by representing the rate of energy 
release per unit area as a peak value over a burn time interval in terms 
of properties Q” and L. However relatively small variations in these 
properties can have a significant effect on fire growth for especially thin 
combustible materials. Overall, the simulation model executed for the 
13 tests yielded fair to good results in predicting flashover. In general, 
the simulation model offers (1) an illustration of how to use fire 
property data for prediction fire growth scenarios, (2) a basis for 
elucidating needed research for improving fire growth models, and (3) a 
preliminary basis for assessing the fire hazard of materials. 
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APPENDIX: PYROLYSIS AREA FORMULAS 


(1) tigo St Sligo + ty: The ignitor burner has caused ignition of region 
(Xp.0, Yp.o) but has not burned out yet. 
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(a) And y,<H, (only the wall is pyrolyzing, Fig. Ala): 
Ap = 2[ yeXp.o + (Xp — Xp.o)Vp.0 + O-5(Yp — Yp.o)(Xp — Xp.o) | 
(b) Or y,=H, (wall and ceiling is pyrolyzing, Fig. Alb): 
Api = 2[F%,.0+ (Op — Xp0)p.0 + 0-5(%p — Xpo)(H — Yoo) 
Zp = Xp — Xp(tu) 


Aen =2| (yp — H)(0-08H) + 325(p— H) 


— 1(0-08H + zy(2—*) |, z,>0 (Ad4) 


P.0 


approx. overlap correction 


OT 


Apy, = 2(yp — H)(0-08H),  z, = 0 
Apc: = min {(2/4)(y, — H)*, wd} 
A, = Api + Apy, + Apcy 


H=2.4m 


Fig. Ala. ¢,,,.=tStig.+t. and y,<H. 
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Fig. Alb. ¢,,,=tStigo+%, and y,=H. 


(2) t>tgot+t. (the ignitor region has burned out and a burn-out 
region (y,, x») Is now computed): 


(a) And y, =A (only the wall is burning, Fig. A2a): 
Ap; = 2[y,0p.0b Apa tee) Yee. t 0 (yaa pclta yt p0)) — (A8) 
Ap = 2[ YoXp.0 ze (x, a Xpo)¥o.0 + 0-5(y. he Yp.o)(Xp ae Xp0)! (A9) 
Ap = Ap) me Ap? (A10) 


(b) Or y,>H,y =H (ceiling is pyrolyzing, but has not burned out 
yet) (see Fig. 4): 


Ap = 2[Hx,.0 = (x, [~ Xp.c)Vn6 + 0-5(H — Yp.o)(Xp vy Xp.o)] (A11) 


Ap? ow 2[ YeXp,0 x (Xp th Xs. OVER + 0-5( ye :* Ypo) (Xb = Kon) | (A12) 
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Fig. A2Za. t>1,..+,. and y, =H. 


Ant 2 0, — H)(0-08H) + 2Zp(Yp — H) 


— 1(0-08H + -(=—*) , z>0 (A13) 


P,o 


ee ee 


approx. overlap correction 
or 
Apy: = 2(y, — H)(0-08H), z,=0 (A14) 
Ap = Ap; — Ap2 + Apy, + Apc (A15) 


(c) »>H, y>H (wall and ceiling have burn-out regions, Fig. 
A2c): 


Ap, = 2[HXx»,.0 ote (xp Hi Xp.o)Yp.0 x5 0-5(x, at Xp,o)(H ot Tel (A16) 
Ap2 = 2[Hx_.0 + (Xb — Xp,0)Vp.0 + O°S(H = Yo0)(Xp —Xpo)] (A17) 
Apy, = by eqns (A13) and (A14) 
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Fig. A2c. t>tigo+ty, and y,>H, yy> H. 
ao 
Zp = Xp — XvG (tu) = Xp — Xo.H (A15) 


Arn =2| (jp H)(0-08H) + 324( 9 ~ H) 


— 3(0-08H + 2) ( —r2) z>0 (A16) 


Xp 
Hf — Yp.0 
Or 
App = 2(y% — H)(0-08H), 2=0 (A17) 
Apcr2 = min {(2/4)[y, — H)? — (¥» — H)’], [wd — 2/4(y,- H)*}} (A 18) 
Ap = Ap; — App + Apn — Apsz + Apciz (A19) 


D.2—Fire Growth Model Source Code 


The FORTRAN source code for Quintiere’s fire growth model is presented below. The 
code used for the analysis in this report is slightly different from the code listed in other reports. 
Changes to the code were made in order to provide a more robust filename management, provide 
more user defined ignition burner output parameters and to increase the accuracy of the main 
program loop. Most of the changes are listed in the source code with the following: [SED] 


Kae Kak kK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK KKK kK Kh hk kkk kK KKK kkk kK KK KKK KKK KKK 


PROGRAM: Fire Growth Model 

VERSION: 1.2 

AUTHOR: Dr. James Quintiere, University of Maryland, 
Department of Fire Protection Engineering 

REVISION: Currently being revised by S. E. Dillon 

DATE: jd — 26-97 


KKK KK KKK KKK KKK KKK KKK KKK KKK KKK KKK kK Kk KKK KKK Kk KKK KKK Kk Kh Kh Kk kkk kk Kk KKK KKK KKK 


OFOFQFOROY OT Ose 


IMPLICIT REAL (a-z) 

INTEGER NMAX, n, burntime, lengthfilename, qburn, qtime, qwidth 

CHARACTER MATNAME*40, name*12 

CHARACTER* 1filename (12) 
c Changed the maximum time to 3600 from 1300. 

PARAMETER (NMAX = 3600) 

PARAMETER (PI = 3.14159, SB = 5.66E-11) 
Cc SB = Stefan-Boltzman constant 

DIMENSION yp(NMAX), yb(NMAX), xp(NMAX), xb(NMAX), t(NMAX), 

* ts(NMAX), time (NMAX), tsO(NMAX), ap(NMAX), qburn(NMAX), 

* qtime (NMAX), qwidth (NMAX) 


VARIABLES 
qig Ignitor energy release rate in kW 
gqigw Ignitor energy per width in kW/m (based on flame height [SED] ) 
qa Total heat release from the pyrolysis area 
q Total heat release {qa*ap + qig} 
yf Flame length 
YP(), Vertical pyrolysis position 
bn Vertical burnout position 
eS: Lateral half-width pyrolysis position 
Cle Lateral half-width burnout position 

Gas temperature 

Surface temperature 
F Surface temperture due to ignition 
), Time from inception of ignition 
ap(), Pyrolysis area 
HCONV(), CONVECTIVE HEAT TRANSFER COEFFICIENT (kW/m**2 K) 
User input data ('filename'.in): 
MATNAME MATERIAL NAME OR ID (40 CHARACTERS MAX.) 
rehlg alee Ambient temperature in Deg. K 
cl Gas parameter 
xrcC Thermal inertia (k*rho*c) 
ho Vent height inm 
w0 vent width inm 


OOOO OLOTOC OO Os, OOF AsO, OSLO O50.) O20 
ct 


621 


OeQeO- O00. 2-50) OG 0 -O40709 sO" CO O80) O70: G)0 


OG sOSO A200) 


Cc 


qf 

qmax 
tsmin 
dtime 
endtime 


Input user determined parameters 
Please note - this input routine is a bare bones, 
method to input data for use while developing the 


Read the 


Height of ignition region inm 

Half-width of ignition region in m 
The number of heat flux levels for the ignitor burner 

The ignitor heat flux levels 


Ignitor energy per width in kW/m: 
The times associated with each ignitor heat flux 


Surface area of roomiin m**2 
Flame net heat flux in kW 


Heat of combustion in kJ/g 


Heat of gasification in kJ/g 
Configuration radiative incident heat flux in kW/m**2 
Tolerance on iterations 


Parameter for ceiling area 
Ignitor incident heat flux in kW/m**2 


Flame length coefficient 


Flame length power 
Total energy per unit area in kJd/m**2 


(ee 


(replaces qig) 
(replaces qigw) 


[SED] 


[SED] 
[SED] 


itqfixig)-—(SB* (tig) A441} 


= 2, 


Ignition temperature in Deg. K 
Lateral flame spread parameter 
Room height inm 


Room width inm 
Room depth in m 


wall; 


r= 4, 


Incident flame heat flux in spread in kW/m**2 
stops run in kw 
Minimum temperature for spread in Deg. K 
Time step in Sec. 


Maximum fire size, 


Maximum run time for the model 


PROGRAM - Fire growth on room surfaces 


a a a a a a ae a a a 


WRITE ~*) 


WRITE 


20 


WRITE *; 


( 
( 
WRITE ( 
( 
( 


Kee 
, 


READ (*, 25) 


oo 


FORMAT (8a1) 


"Enter 
(filename (i), 


algorythms. 


"filename": ' 


i=1,8) 


"filename'.IN to be used 


Determine the length of 'filename' 
DO 30 I=8,1,-1 


LE 
a 


30 CONTINUE 
35 lengthfilename = I 
Open the input file 
CALL FILE(filename, lengthfilename, '.in',3,name) 
(1, FILE=name, STATUS='OLD', ERR=2000) 
Read the data from the input file. 


OPEN 


READ (1, 
READ (1, 


50) 


si 


(filename (I) .NE.' 
(1. EO. A+) 


GOTO 5000 


MATNAME 


ting 


” 


GOTO 35 


[SED]. 


622 


(sec) 


[SED] . 


[SED] . 


[SED] 


ceiling) 


[SED] 


quick and dirty 


READ Ces Cor 
READ({L)) y= xre 
READ (CL; == ) “He 
READ(1, *) w0 
READ(1, “)syveo 
READ (i; “*“}"xpv 


READ(1, *) qnum 
DO. (45enc= yi) qnum 
READ(1,40) qburn(n),qwidth(n),qtime(n) 


FORMAT (317) 
CONTINUE 

READCUL . ees 
READ(1, *) qfnet 
READ GIS? *) Fhe 
READ Ae * ) exe 
READ(i,; *)eqrizg 
READ(L;. -*) stor 
READ (= )t 
READ(, *)earexig 
READ(i;7 7) c2 
READ (i; = 0CS 
READ(1, *)c4 
READ (1) **) “tig 
READ(1, *) phi 
READ (Lp *) th 
READ (lL, *) w 
READ (i, °*) +d 
READ(1, *) q£ 
READ(1, *) qmax 
READ(1, *) tsmin 
READ(1, *) dtime 
READ(1, *) endtime 
CLOSE (1) 


Initialized variables: 

Time from aplication of ignitor 

time(1l) = 0 

Room gas temperature in Deg. K 

t(1) =F293 76 

Avg. room surface temperature in Deg. K 

ts (1)?= 729330 

Surface temperature due to ignitor in Deg. K. 
ts0 (1) ="298.0 

Pyrolysis area in m**2 


ap(1) = 0.0 

Pyrolysis front vertical coordinate inm 
yp(1) = 0.0 

Burnout vertical coordinate inm 

yb(1) = 0.0 

Lateral pyrolysis coordinate inm 

xp(1) = 0.0 

lateral burnout coordinate inm 

xb(i) = 0.0 


The ignitor burner heat flux level counter (1 = initial) 
burntime = 1 

Energy release rate in kW 

gig = qburn(burntime) 
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Q 


gigw = qwidth(burntime) 
q = qig 

Index 

kount = 0 

Index 

kyp = 0 

kyb = 0 

Ceiling jet dummy values 


xph 
xbh 


= 0. 
= 0. 


Ignitor ignition time with no heating in Sec. 
tmigO = 1.E06 


Open the output file 

CALL FILE(filename, lengthfilename, '.out',4,name) 
OPEN (1, FILE=name, STATUS='UNKNOWN' ) 

WRITE(1, 50) MATNAME 


* 


* 


WRITE (1, *) "time(s) yp(m)  xp(m) yb (m) xb (m) apm 
* (mac) t (K) ts (K) q( kw)" 
WRITE (1,ya000)) time (th mye tl) Fixp(l),cybil); xbily, papi{l), 
teippets(1l), gq 
50 FORMAT (A40) 
MAIN LOOP 
Advance Time 
DO 100 n = 1, endtime 
Increase time step 
time(n + 1) = time(n) + dtime 


Does the ignition burner need to be adjusted? [SED]. 


IF (time(n + 1) .EQ. gtime(burntime)) THEN 
Ignition burner is adjusted based on the input file [SED] 
burntime = burntime + 1 
qig = gqburn(burntime) 
gqigw = qwidth (burntime) 

Has the maximum burner time been reached? 
IF (burntime .GT. qnum) THEN 
Ignition fire is turned off 
gig = 0. 
gigw = 0. 
END IF 
END IF 


Does Ignition Occur? 


ie 


(time (meen oenLT.. vemiig0) aaTHEN 
Set burn region to zero 
yp (rit 1) se 
xp(n +) == 901 
yb (fiat 1) =O 
xb (n+) eau 
ap(n (aL. 230% 
Compute room gas temperature due to ignitor 
CALDGtempitinet 1), t(n), ap(n + 1)q) *#imeim +#01)3 
lic, #xliginet p qrfig,,.qigyatintsew0j;1h0, as, 
ere ely, tol) 
Compute room surface temperature 
CALL Stemp(ts,t,time,dtime,0.,0.01,xrce,tinf,tol,n, NMAX) 
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a Compute surface temperature at ignitor 
CALL Stemp(ts0,t,time, dtime, qflxig,0.0,xrc, tinf, tol,n, NMAX) 


(e Does ignition Occur? 
EF. (tsOin +*:l); .GEs> t1g) THEN 
c Record ignition time. 
tmigO = time(n + 1) 
c Compute the energy release rate per unit area (kW/m%2) 
qa = (he/xl) * (qfnet + qrfig + sb * t tuieal) seep 
e Compute initial burnout time. 
tmbO = c4/qa 
c Set initial values of pyrolysis region 
yp(n + 1) = yp0 
xpin + 1) = xpd 
Vii Tiers: ae) eee 
p39" Cancer a at ae 6 Re 
apne Teac. | A xXDO UAT VypU 
ELSE 
GOTO 410 
END IF 
ELSE 
200 CONTINUE 
c Controis burn calculations 
cl) . Tee? 2 ee eee eee 


IF (time (n + 1) eGT. (Gomig0) + stmb0)) eS THeN 
250 CONTINUE 
S Ignitor burnout has occurred. 
fe This section computes burnout fronts (yb & xb) as 
C well as pyrolysis fronts (yp &S@xp) forgntiiby 
S Runge Kutta 
G 
C 


Séteinitial conditdonssforgaeyb andexb 
kount = kount + 1 
Lee kount sr. BOs. )) THEN 


yb(n) = yp0 
xb(n) = xp0 
END IF 
260 CONTINUE 
(e Calculate YPE 
CALI CRT £11, .Vpi( nips VOD) aint (je ms (ia), 
ne, xl, qinet, qriig, C2, C37 adiows; ero; Lao 
ype = yp(n) + £11 * dtime 
c Is the surface temperature less than what is needed 
S for lateral flame spread. 
IF (ts(n) .LE. tsmin) THEN 
= No lateral flame spread 
xpe = xp0 
xbe = xp0 
ELSE 
c Calculate lateral flame spread 
CALLI FSX ESigmes (nn), 
cs crc, tip; aphay 


XpeV= "XD (nN) ors.) eC ime 
CALL .F4(£41,.xpi(n), xiao acim) > 
a he,’ xl,oginet,) qriiggi cot) 
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xbe = xb(n) + £41 * dtime 


END IF 
= Calculate ybe and ape 
~. CALL F2 (E21 yp ieee on) Fo) 7 
i he, Vk), Sqinet) “orlig?? cf) 


ybe = yb(n) + £21 * dtime 
CALL Area(ape, ype, ybe, xpe, xbe, yp0, xp0, h, w, d, 
S Tecmo eeu, ctame(n + 1) ,\4xpn/iixbh) 
CALL Gtemp((te7sti(n)s ape, time(n +41), hc, -xl, -qfnet, 
= (mei, dig A-ctink, wo, nO as,-xre,; 'cl1, tol) 
e Estimated gas temperature in iteration 
t(nat 1) ee2 ce 
CALLYStemp (ts, t, time, dtime; 0". ,0.01, xrc;tint, tol,n, 
= NMAX ) 
c Estimated surface temperature in iteration 
tse = ts(n + 1) 
c Complete Runga-Kutta for yp, yb, xp and xb 
CALLISi (fiz, ypey*ybe;}) teyrtse, 
. Deer pect nee pe oriiGg,@C2z,1 C3, digw, “xrc, tio}. qty 
Vptnise) = yo(n) + (fla aie) Serduame 97'2. 
CALL F3(£32, tse, 
*¢ xeC, Eig, phi) 
<Dimereely= xp (in) + (£31 4+ £32) * dtime / 2. 
CALL F2(£22, ype, ybe, te, 
Rojee, soEnet, Grtig, ¢cé4) 
omer eaye=nybin) ye hi2l, +8 £22) )*, dtame* /) 2. 
CALL F4(f£42, xpe, xbe, te, 


= Ro,exl, ginet, «ritig, c4) 
Soe) = xb (n) + (£41 + £42) * dtime / 2. 
ELSE 

300 CONTINUE 
S Ignitor burnout has not occurred 
e This section computes yp and xp 
c by a second order Runge Kutta method 
e Set burnout region to zero 


VoAT 1) = QO. 
eb a+) 1) Sag 


C Compute estimated (ype) 
CaGorei(ril, yo(n); yo(n), tin), tsi{n), 
= Noe, ext, .aqinet, qriigq; ce, \e27e-qigw,- xnc,? tig?! qf) 
ype = yp(n) + £11 * dtime 
o Compute estimated (xpe) or fixed xp 
Lie its (ny. eLeSetsmin). OTHEN 
c No lateral flame spread 
xpe = xp0 
ELSE 
re Calculate xpe 


Chibereotgs., tsin);, xxrc,-tig,. phi) 
xpe = xp(n) + £31 * dtime 


END IF 
: Compute estimated pyrolysis area (aps) 
CALL Area(ape, ype, yb(n + 1), xpe, xb(n + 1), yp0, xp0, 
* hptw, ed; er, etmigd,-tmb0,;’time(n’+ 1); xph, xbh) 
cs Compute estimated gas temperature (te) 
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400 


CALL Gtemp(te,,t(n), ape, -taime(m + 1) » hes) lr ginet, 
qrfig;. gig,#tin£;,,w0,7h0, nas,*xrc, ecw etc 
Compute estimated surface temperature (tse) 


tlpete ) ase 
CALL Stemp(ts,t,time,dtime,0.,0.01,xrc,tinf,tol,n, 
NMAX ) 


tse = ts(n + 1) 
Compute yp(n + 1) 
CALL WEL (£12, ype, yo(netelds> ters tse, 
he »exd,,wqinet, axnfignn(C2Z,.163) gidwyn xrc, tigs ate 
yp(ni t+: 1) .= ye(ndetta(thy t9f12)) peedtime: /ei2s 
compute xp(n + 1) 
CALL, -F3%£32,.<¢se, 
¥rC, BCi0; en.) 
xp(n + 1) = sep (meet £31 ehpfe2lecGedt imevier2. 
END IF 


CONTINUE 
Check if yp and yb have reached ceiling 
Compute values xph and xbh 
IF (yp(ntl1) .GE. h) THEN 
kyp = kyp +1 
END IF 
IF (kyp .EQ. 1) THEN 
xph = xp(n+1) 
END IF 
TF (yb (n+Lje).GE. Bb) TEEN 
kyb =" kybDetae 
END IF 
TE o(kyb. (Sh Os asl) ee HEN 
xbh = xb(nt+1) 
END IF 

Computewap, ity d, «tS, at Benepe) 

CALLyvArea(ap(n +t 1); ypiln ey lis vDAnetiaexpineteae), 
xbin + i), vpQjampoy dypewi eee) omo0 atmo; 
me. ( tal) > xpi, exon) 

CALIn Gtemp(t(n + 1), t(n), ap (nosed jos timeims te cane, a, 

qfnet, qrfig, 

dig; stint, .w0,, Dog as, xrCy aC. sero 
CALE.Stemp(ts, t, time, dtime, 0.,e020l @re,cint; tolgn: 

NMAX ) 

[Pints (ne+ fae . beet smi pet EN 

No lateral flame spread 
IF (time(n. + 1) ».GEee (tmigO “.tmb0))) THEN 
Burnout occurs 
xp(n + 1) = xp0d 
xb (ni +h) ¢aexp0 
ELSE 
Burnout does not occur 
sO + 1) = Exp 
xb.(nearbalhy 240% 
END IF 
Recompute ap, t and ts 
CALIGATrea (apn, thd )iyacyp (Nites! )wi e¥D (natty L)eoxpinet LA 
xbilnechs des -yp0juxp0 O.hyew,) ad, ir, emig0, tmb0, 
time(n +1). oph7s xbh)} 
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CARLAGi empty net ele ene seapine +2), time(n +1). he, 


xly Ginet;- criig, 
Gigemusnt, wo, hO, as, xre, cl, tol) 
CALE JStemot tat, 41m], ceime, 0O., ,0.01,;xre,tinf;tol,n, 


NMAX ) 
END IF 
C Compute heat release 
410 CONTINUE 
e Compute the energy release rate per unit area (kW/m%2) 
gqai=e(hc ¥ S0)* *miotnetesiqrnéig: sesh, siti(nidte 1) +44) 
c Compute the total heat release 
gq = qas* ap(n +°1) + qig 
c Has the maximum fire size been reached? 
IF ( q .GT. gqmax) THEN 
c End the program 
write(1, 1000) time(n+1l), yp(n+1), xp(nt1l), yb(nt1), 
xDinviyyeapunia);)t intl), ts(ntl), g 
CLOSE (1) 
WRITE (*,*)name,' created.' 
GOTO 20 
STOP 
END IF 
END IF 
Write (ape. OO0Nm time (n+), ypi(ntl); xp(nt1l), ybo(ntl), 
cet apn nt 1). tuntljgets (nt1) og 
S Continue time advance 
100 CONTINUE 
500 CONTINUE 
CLOSE (1) 
WRITE (*,*)name,' created.' 
GOTO 20 
STOP 
ere ORMAT (Fi 27 okee On cgioeky FO. 3, 3X, F5.3;73X,F5.3, 
Mens So Shee) 2, 3X, E722, Shp EI st) 
2000 WRITE(*,*) 'ERROR -- Could not open ',name 
GOTO 20 
5000 STOP 
END 
c SUBROUTINES 
So SS MED Gen par een 
SUBROUTINE FILE(filename, lengthfilename, extension, lengthext, name) 
c [SED] 
c This subroutine is used to create the appropriate filenames that 
c can be used by the main program for opening the .in and.out files. 


dimension filename (8) 

dimension extension (4) 

character filename, extension, name* (lengthfilename+lengthext) 
length = lengthfilenamet+tlengthext 


Fe 

DO 5 i=lengthfilename+1, length 
filename(i) = extension(j) 

j = jtl 
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5 CONTINUE 

WRITE (name, 10) (filename (I),I=1,length) 
10 FORMAT (12al) 

RETURN- 

END 


SUBROUTINE, FLI(fOL, syp, avd; t, ts, nc, xl, eoinet,.artsd,ec. ce, 
bs gigw, xEc,. big;, of} 
Computes RHS function of first order DE for yp 
PARAMETER (SB = 5.66E-11, PI = 3.14159) 
Energy release rate per unit area (kW/m**2) 
ga = (hew/. xlje *edqineta+ qriig +eShace2 25 
Flame length (m) 
IF ( yb ,LE. c2i* gigw**c3ja THEN 

Vies=e Voet C2 * {(qigqwit'qa 47 typos yb)j¢ see 
ELSE 

yi = yb +i clm*eiqatt*h vpn —-. yo)oFGeC3 
END IF 
Ignition, times (Sec: } 
tmigue= (Piv/.4 je *exrcets (tig =) ts) * Zen @ciea 
FOlLL= (vi. - yp) /, tmig 
RETURN 
END 


SUBROUTINE. F2(f02, yp, yb, t, Rc; xiGp cine, soriig,ecs 
Computes RHS function of first order DE for yb 
PARAMETER (SB = 5.66E-11) 

ety Es Whale iil Salyw tl ilepae(che. Co lebahaellop im, ys; (A (etter) 

tmb = c4 / qa 

pO VAT EST Tek sy eg Gi: Seale, 

RETURN 

END 


SUBROUTINE F3(£03, ts; (xrc,, tigpepn.) 

Computes RHS function of first onder DE for xp 
£037 = phi if [xrc ese (tidgee mts) =*2) 

RETURN 

END 


SUBROUTINE. F4(.f047.xp,) xb,ant, gicwesi, cinet,; .orrig;, «c4) 
Computes RHS function of first order DE for xb 
PARAMETER (SB = 5.66E-11) 

qa = (he / x1) .% slatnet atworiiour SB: =r 4) 

tmb = c4 / ga 

£04 = (xp - xb) / tmb 

RETURN 

END 


SUBROUTINE, Area (ap,..yp,- YDr4 XP, xb, yp0, xp0, h,' w, d, r, 
- tmigO, tmb0O, time, xph, xbh) 

Calculates the fire growth on a wall and ceiling 

Wall pyrolysis area with no burnout, apl 
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Wall burnout area, ap2 
PARAMETER (PI = 3.14159) 
timeb = tmigO + tmb0 
IF (time .LE. timeb) THEN 
IF (ypv.LEs h} \THEN 
No burnout and only wall pyrolyzing 
apugelee Cfo eeer xp —-xp0) * -yp0 F#€0.6f4.49p =19p0) 
“ * I AKXDU=" xp Uy) 
ap = apl 
ELSE 
No burnout and ceiling is pyrolyzing 
abe cee ee ecco (xp - xp0) *-yp0 i+ (OWS e7 (hvexyse) * 
- (xe ox pUy.) 
Wall ceiling jet 
Zp: =ExD) -exph 
IF (zp .GT. 0.) THEN 
Abjeeoec mara Yousenio* 0. 06urt ho+e0r>). teegp: > 


’ (yp) -) hh) =@0e 5e%.GOt08*h #0zp) ***2 * 
* (xp-xp0) / (h-yp0) ) 
ELSE 
Be ieemesw ts RVD." fh), * 0.08 * nm 
END IF 


Ceiling area 
Boece r. 7 £) * (yp "=8h)}*A2 
aclmax = w*d 
apcl = AMINI1( apcl, aclmax ) 
ap = apl + apjl + apcl 
END IF 
ELSE 
LE (9p Lege hy. THEN 
Burnout and only wall is pyrolyzing 
acter oeee ca type tnxpOl? Yxpieexp0) *)3p0 + 0.5 * (yp - yp0) 
4 eo = xp0)') 
Beemer (YD 8 xo0 + (xb = xp0} * yoO + 0.5 *) {yb -— yp0) 
* (sO @= x0 ),) 
ap = apl -‘ap2 
ELSE 
IF (yb .LE. h) THEN 
Wall burnout and ceiling is pyrolyzing 
Apomer 2.  eih) * sposas(xp <expd)) Beypo oO 752% 


= A> pOoyier? .xpe -exp0)) 7 
Spomeme ss (vbr 2txp0 ie txb — xp0) * ypO + O85 "* 
i oem pOjae Atxbesixps= )t) 


Wall ceiling jet 

Zp, = xp - ph 

LP etzpe GT.) 0.) . THEN 

epyigeez.) * (yp born) 20.108 Bshpes0"gs * zp.* 


* (yp - h) - 0.5 * (0.08*h + zp) **2 * 
3 (xp-xp0) / (h-yp0) ) 

ELSE 

api le= 2709 wp = Sh) e* “0.308 Ah 

END IF 


Ceiling area 

apCle = to ate vp em hh) **2 
aclmax = w*d 

apcl = AMIN1( apcl, aclmax ) 
ap = api — ap2 + apjl + apcl 
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20 


ELSE | 
Wall and ceiling have burnout portions 
apl = 2. * (h * xp0 +. (xp -— xp0) .*);yp0 + 0Geme SAGs ia0y : 


* Fae Ck De XD) 
ap2.= 2. * ({h * xp0.+ (xbi- xp0) * vp0 7GRS te (Gee vee 
* Re (xia DO) 


Wall ceiling jet : 
20a ee 

TE (2p gaCTe Oo THEN 

2p3 de 2 et br i Oe eh area Ze 


> (yp — hh} oen0.9) * e008 Fhutmanis 2 2ae 
4 (xp-xp0) / (h-yp0) ) 
ELSE 
AO eee Oe en Vek) i) oot eam 
END IF 


zb = xb - xbh 
LF (2b oGT.@@.) THEN 


apj2.=.2.* ((yb.~ Bys*On08. 4*50100.5 Bieb=* 
‘ (yb =n) 
* - 0.5 * (0.08 * h + zb)**2 * (xb - xp0)/ (h-yp0)) 
ELSE 
apj2 = 2. * (yb - h) * 0.08 * h 
END IF 


aclZmx = 'w * d= (Pi/r) * (yee) 22 
ape le = Hele Ge Vpn eee ae yan ear) 
apcl2 = AMINI(apcl2 , acl2mx) 
ap.= ap! —~yapzet api l= apy Zatmapcre2 
END IF 
PND mess 

END IF 

RETURN 

END 


SUBROUTINE Gtemp(tnl, t, ap, tame, he, xl; qinety Gauri; 
Gig acini, wopen .1eS, ok LCheC i neo) 
Computes room gas temperature from correlation by iteration 
PARAMETER (SB = 5.66E-11) 
Cerne ae 
First new value of t by Gauss Siedel iteration 
ga = (he / x1} * (qinet + qrfiig tiese eet 44) 
G= Gigs* caeegau 
Modified to account for corner firejicoe ties. Gs foe a 
P2beetant, * (te Ot 2a al (c/s lc ee es 
PO FL Se Se) ie ae meee Chie cd ck a 


* Las) te SORT (SEC! Ja Lame) ji) eet a) mero ee 
estimate new value of t by Regula Falsi method 
CONTINUE 

gqal.= (ne l/ xl).*% (atnet + oqrfig,+ Shee tise 


Gle= gig +) cad * ap 
fl = tint * (4.0 + 2.2 -* (qh belo tit 0 


a hO**T 5) YF 26. /6 Saas al Clie wo ene aS eee 
* (as. * SORT ixrco'/ time) ))**:(1. 773-0) )) eet 
new tl 

rele ay 


qaZ_= (he / xl) * (qtnet_4 oriigy tase om t2 ee) 
92 = (01o Tt eGazen ap 
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Pie CI ee? ee ee Le ae ae 
DO Wat kemce Sater cL WO TOS eo) / 
F (as * -SORT(xrc lyme os / 330) )) = 62 
S new t2 
E26 tlertie(tZee bl) PRABS( Di velf ABS( ft), 2 ABS(£2)»,) 
Drei ABSit? - tl) /°t2 2GE. tol) THEN 
GOTO 20 
END. IF 
tnl = t2 
RETURN 
END 


SUBROULINEMSt embiwo eG come, dtime,; “qfixirg;"HCONV,xrc;tint,; 


= tol, n, NMAX) 
c Computes average room surface temperature by semi-infinite wall 
c model solving integral equation at surface by iteration 
é This version includes convection 
PARAMETER (SB = 5.66E-11, PI = 3.14159) 
DIMENSION ts (NMAX), t(NMAX), time (NMAX) 
wi = eS (hy) 
sum = 0 
re {my .@q. 2) THEN 
sum = 0.0 
ELSE 
DOS ies Ll, n= 1 
pepemufaxagi+ sb* (ti(5g)**4 - ts (p)a*4) 
* Sereonv=(t(4)—ts(j))) / 
= Somer came(n + 1) — tame(4) ) 
Pomerat Sb * (t(j oF 1) **4 = ts (4+ Ty *4) 
* BONS TAL ee rset) ) 
* See tame (n +1) — time (7 + 1)) 
Sone esumet £1 + £2 
100 CONTINUE 
END IF 
Te Pe te ete She (tin) #44; —-tsin)**4) 
x te BION Veit (xjde= it Sm)y) 
10 CONTINUE 
feces sgl lor oats SBio* atin Gatelatt em atlas 4) 
i hRBCONVA (tint lent) 
Ceeehtins  tatcum dt imate 2 o Out £61 oto £62 jour 
* Semi Gatime )jn/ SOrt(Pl ter) 
TE, WARS Atl ie/at2 .sGE. tol) THEN 
tl = t2 
GOTO 10 
END IF 
tSiot i st2 
RETURN 
END 
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